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Abstract 

We present a lattice calculation of the electromagnetic form factor of the pion 
obtained using the tree-level Symanzik improved gauge action with two fla- 
vors of dynamical twisted Wilson quarks. The simulated pion masses range 
approximately from 260 to 580 MeV and the lattice box sizes are chosen in 
order to guarantee that M-^L > 4. Accurate results for the form factor are 
obtained using all-to-all quark propagators evaluated by a stochastic pro- 
cedure. The momentum dependence of the pion form factor is investigated 
up to values of the squared four-momentum transfer ~ 0.8 GeV^ and, 
thanks to the use of twisted boundary conditions, down to ~ 0.05 GeV^. 
Volume and discretization effects on the form factor appear to be within 
the statistical errors. Our results for the pion mass, decay constant and 
form factor are analyzed using (continuum) Chiral Perturbation Theory at 



next-to-next-to-leading order. The extrapolated value of the pion charge 
radius is (r^)^^^^ = 0.456 ± O.OSOgtat. ± 0.024 

syst. ™ '^i^'^ agreement with 
the experimental result. The extrapolated values of the pion form factor 
agree very well with the experimental data up to 2± 0.8 GeV^ within 
uncertainties which become competitive with the experimental errors for 
Q2 > 0.3 GcV^ The relevant low -energy constants appearing in the chiral 
expansion of the pion form factor are extracted from our lattice data, which 
come essentially from a single lattice spacing, adding the experimental value 
of the pion scalar radius in the fitting procedure. Our findings are in nice 
agreement with the available results of ChPT analyses of tt — tt scattering 
data as well as with other analyses of our collaboration. 



1 Introduction 



The investigation of the physical properties of the pion, which is the hghtest bound 
state in Quantum Chromodynamics (QCD), can provide crucial information on the 
way low-energy dynamics is governed by the quark and gluon degrees of freedom. 
In this respect for space-like values of the squared four-momentum transfer, = 
—<f > 0, the electromagnetic (e.m.) form factor of the pion, F^.^Q'^), provides 
important insights on the distribution of its charged constituents, namely valence 
and sea light quarks. At momentum transfer below the scale of chiral symmetry 
breaking {Q^ < 1 GeV^) the pion form factor represents therefore an important 
test of non-perturbative QCD. 

The current experimental situation is as follows. For values of < 0.2 GeV^ 
the pion form factor has been determined quite precisely at CERN SPS [I] by 
measuring directly the scattering of high-energy pions off atomic electrons in a 
fixed target. At higher values of the pion form factor is extracted from cross 
section measurements of the reaction ^H{e, e'7i~^)n, that is from electron quasi- 
elastic scattering off virtual pions in a proton. The separation of the longitudinal 
and transverse response functions as well as the extrapolation of the observed 
scattering from virtual pions to the one corresponding to on-shell pions have to 
be carefully considered for estimating the systematic uncertainties. Using the 
electroproduction technique the pion form factor has been determined for values 
in the range 0.4 9.8 GeV^ at CEA/Cornell p], for = 0.35 and 0.70 GeV^ at 
DESY [3], Hj and, more recently, for in the range 0.6 1.6 GeV^ [5j and for 

= 1.60,2.15 and 2.45 GeV^ |6] at the Thomas Jefferson National Acceleration 
Facility (JLab). A careful reanalysis of the systematic uncertainties for the data 
of Refs. [31 m E] has been carried out in Refs. [3, IB]- 

It is well known that at small values of the pion form factor can be repro- 
duced qualitatively by a simple monopole ansatz inspired by the Vector Meson 
Dominance (VMD) model with the contribution from the lightest vector meson 
(Mp ^ 0.77 GeV) only. This is not too surprising in view of the fact that in the 
time-like region the pion form factor is dominated by the p-meson resonance [S]. 

More interesting is the quark mass dependence of the pion form factor, which 
can be addressed by QCD simulations on the lattice and by Chiral Perturbation 
Theory (ChPT). The latter, which is known at next-to-leading (NLO) [10] and 
next-to-next-to-leading (NNLO) order for the pion form factor, can be used 
as a guide to extrapolate the lattice results from the simulated pion masses down 
to the physical point, obtaining at the same time an estimate of the relevant low- 
energy constants (LEC's) of the effective theory. 

Initial studies of the pion form factor using lattice QCD dates back to the 
late 80 's [121 [13] giving strong support to the vector-meson dominance hypothesis 
at low Q^. Within the quenched approximation, which neglects the effects of 
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the sea quarks, several lattice investigations have been carried out using Wilson 
[H], Sheikholeslami-Wohlert [15], twisted Wilson |T6J and Ginsparg- Wilson [17] 
fermions. The effects of the quenched approximation might be limited because, 
thanks to charge-conjugation and isospin symmetries, the e.m. pion form factor 
receives no contribution from the so-called disconnected diagrams in which the 
vector current is attached directly to a non- valence quark (see Ref. [IS])- However 
there are effects from sea quarks which do not interact directly with the external 
current, and they can be taken into account only by performing unquenched gauge 
simulations. 

There are few results for two flavors of dynamical fermions from JLQCD [TH] 
and QCDSF/UKQCD [11] collaborations adopting Clover fermions and again from 
JLQCD [20j using overlap quarks. Finally only two studies with three flavors of 
dynamical quarks are available to date, namely from Ref. [14j, where domain- wall 
valence quarks and Asqtad sea quarks are mixed, and from Ref. [21], where the 
domain-wall formulation is used for both sea and valence quarks. 

As far as the lattice results for the (squared) pion charge radius at the physical 
point are concerned, the present situation is a bit puzzling. Some collaborations 
[m [18] have found that their extrapolations underestimate significantly (up to 
~ 30%) the well-known experimental value (j'^y^P- = 0.452 ib 0.011 fm^ [22], while 
other collaborations [121 [201 121] have obtained values consistent with experiment 
within the errors. 

The European Twisted Mass (ETM) collaboration has recently produced a 
large number of gauge configurations with two flavors of dynamical quarks [2^1 [211 
[25j using the Wilson twisted-mass fermionic action [26j and the tree-level Symanzik 
improved (tlSym) gauge action [27J. In order to obtain (almost) automatic 0{a) 
improvement the Wilson twisted-mass fermions have been tuned to maximal twists 
[28] . An intensive, systematic program of calculations of three-point correlation 
functions relevant for the determination of meson form factors both in the light 
and in the heavy sectors has then been started. Preliminary results, concerning 
the vector and scalar form factors of the pion, the Isgur-Wise universal function 
and the transition form factors relevant in Ki^ and D 7t{K) semileptonic decays 
have been presented in Ref. [2P] . 

In this paper we concentrate on the e.m. form factor of the pion and we present 
the results of several measurements performed with pion masses in the range from 
~ 260 MeV to ~ 580 MeV, using six values of the quark mass at a lattice spacing 
of ~ 0.09 fm and two values of the quark mass at a lattice spacing of ~ 0.07 fm. 
The lattice box sizes are chosen in order to guarantee that M^r-L > 4 for minimizing 
as much as possible finite size effects. Thanks to the use of all-to-all propagators 
evaluated by the stochastic procedure of Ref. |^ (see also [24]) the statistical pre- 
cision of the extracted form factor is quite impressive. The momentum dependence 



4 



of the pion form factor is investigated up to values of the squared four-momentum 
transfer ~ 0.8 GeV^ and, thanks to the use of twisted boundary conditions 
(EC's) [an [32], down to ~ 0.05 GeVl The Q^.ghape at the simulated pion 
masses is well reproduced by a single monopole ansatz with a pole mass lighter 
by 10% ^ 15% than the lightest vector-meson mass. Volume and discretization 
effects on the form factor are estimated using few simulations at different volumes 
and lattice spacings, and they turn out to be within the statistical errors. 

The extrapolation of our results for the pion mass, decay constant and form 
factor to the physical point is carried out using (continuum) ChPT at NNLO [11]. 
The extrapolated value of the (squared) pion charge radius is {r'^)P^y^ = 0.456 ± 
O.OSOg^g^^ ± 0.024gyg|^ in nice agreement with the experimental result (r'^y^P- = 

0.452 ± 0.011 fm^ [22]. The extrapolated values of the pion form factor agree 
very well with the experimental data up to ~ 0.8 GeV^ within uncertainties 
which become competitive with the experimental errors for Q2 > 0.3 GeVl The 
relevant low-energy constants (LEC's) appearing in the chiral expansion of the pion 
form factor are extracted from our lattice data adding in the fitting procedure the 
experimental value of the pion scalar radius [TT]. The latter helps constraining 
one the LEC's (£4), which in turn is beneficial for reducing the uncertainties of the 
other LEC's. We get: ii = -0.4±1.3±0.6, h = 4.3±1.1±0.4, = 3.2±0.8±0.2, 
£4 = 4.4 ± 0.2 ± 0.1, 4 = 14.9 ± 1.2 ± 0.7, where the first error is statistical and 
the second one systematic. Our findings are in nice agreement with the results of 
ChPT analyses of vr — vr scattering data [33]. The values found for £3 and £4 agree 
very well both with earlier ETMC results from Refs. [231 [2^ and with the recent 
ETMC determination of Ref . [3l] . This is quite reassuring because different kinds 
of systematic uncertainties may affect the two analyses: the present one being a 
NNLO analysis limited mainly to data from a single lattice spacing, and that of 
Ref. [SI] having two values of the lattice spacing but limited mainly to NLO in 
ChPT. 

The plan of the paper is as follows. In the next Section we briefly discuss the 
implementation of twisted EC's for the quark fields. In Section [3] we present the 
calculation of two- and three-point correlation functions performed in the Ereit 
reference frame, where the values of four-momentum transfer are independent 
of the simulated pion mass. We also briefly show the stochastic procedure used 
for our unbiased estimate of the all-to-all propagators employed in this work. 

In Section m we firstly illustrate the very precise results obtained for the renor- 
malization constant of the vector current and then we compare them with other 
determinations. Our accurate results for the momentum dependence of the pion 
form factor for the various simulated pion masses are presented and both volume 
and discretization effects are investigated. 

In Section |5l using a single monopole ansatz to fit the momentum dependence of 
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the form factor, the charge radius and the curvature are calculated at the simulated 
pion masses and analyzed both in terms of the ChPT expansion at NNLO from 
Ref. pT] and adopting a simple polynomial fit. 

The mass and momentum dependencies of our lattice points for the pion form 
factor are analyzed in Section without any model assumption, but using only 
the functional forms dictated by ChPT at NNLO. 

In Section[7]the final values of the relevant LEC's, including the estimate of the 
systematic errors, are presented and it is shown that the extrapolated form factor 
at the physical point agrees very well with the experimental data in the whole 
range of values of considered. Finally Section [H] is devoted to our conclusions. 

2 Lattice all-to-all quark propagators with twisted 
boundary conditions 

In lattice QCD simulations the spatial components of the hadronic momenta "Pj 
(j = 1,2,3) are quantized. The specific quantized values depend on the choice 
of the BC's applied to the quark fields. The most common choice is the use of 
periodic BC's in the spatial directions 

i){x + tjL) = ^{x) , (1) 

that leads to 

where the n/s are integer numbers. Thus the smallest non- vanishing value oi pj 
is given by 27i/L, which depends on the spatial size of the (cubic) lattice (V = 
L^). For instance a current available lattice may have L = 32 a, where a is the 
lattice spacing, and ^ 2.5 GeV leading to 27r/L ~ 0.5 GeV. Such a value 
may represent a strong limitation of the kinematical regions accessible for the 
investigation of momentum dependent quantities, like e.g. form factors. 
In Ref. [31j it was proposed to use twisted BC's for the quark fields 

ip{x + ijL) = e^^'^^ i^ix) (3) 

which allows to shift the quantized values of pj by an arbitrary amount equal to 
2TT9j/L, namely 

^2n 2n ^2% . ^. 

Pj = Pj + dj^ = ^jj^ + ' 

The twisted BC's ([3]) can be shown [31j to be equivalent to the introduction of 
a U(l) background gauge field coupled to the baryon number and applied to quark 
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fields satisfying usual periodic BC's (the Aharonov-Bohm effect). In Ref. [32] the 
twisted BC's were firstly implemented in a lattice QCD simulation of two-point 
correlation functions of pseudo-scalar mesons. The energy-momentum dispersion 
relation was checked confirming that the momentum shift 2TT9j/L is a physical 
one. In Ref. [33] the twisted BC's were firstly applied to the calculation of the 
vector and scalar form factors relevant to the K ^ it semileptonic decay. It was 
shown that the momentum shift produced by the twisted BC's does not introduce 
any additional noise and easily allows to determine the form factors with good 
accuracy at quite small values of Q^, which are not accessible when periodic BC's 
are considerecll]. 

On the lattice, for a given flavor, the all-to-all quark propagator S{x, y) = 
{ip{x) ip{y)), where (...) indicates the average over gauge field configurations 
weighted by the lattice QCD action, satisfies the following equation 

^D{x,z) S{z,y) = 6^^y (5) 

where D{x, z) is the Dirac operator whose explicit form depends on the choice of 
the lattice QCD actioiJl. In what follows we work with the fermionic twisted-mass 
Lattice QCD (tmLQCD) action with two flavors of mass-degenerate quarks given 
in Ref. [23], tuned at maximal tmst in the way described in full details in Ref. |2^ . 
Therefore, in the so-called physical basis the operator D{x,z) is given explicitly 
by 

D{x, z) = K{x, z) + ^75X3 W{x, z) + am 5x,z , (6) 
K{x,z) = -"^-f^, {6^,z^af, U^{x) - Ul{z)} , (7) 



2 



4 

r 



W{x, z) = (4r + amcrit) 5x,2 " 2 XI {^^•.^-«A U^{x) + Ul{z)} , (8) 

where U^(x) is the gauge link, m is the bare twisted quark mass, merit is the critical 
value of the untwisted quark mass (needed to achieve maximal twist), Ts is the 
third Pauli matrix acting in flavor space, and r is the Wilson parameter, which is 
set to r = 1 in our simulations. 

We now want to consider the case in which a valence quark field satisfies the 
twisted BC's ([3]) in the spatial directions and is anti-periodic in time. This is at 
variance with what has been done in the production of the ETMC gauge configura- 
tions, which include two sea quarks with periodic BC's in space and anti-periodic 



"'^We mention that a new application of twisted BC's to the evaluation of the vector form 
factor at zero-momentum transfer has been proposed in Ref. |21j . 
^We omit in this Section color and Dirac indices for simplicity. 
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ones in time [2^. However it has been recently shown [36] that for many physical 
quantities, which do not involve final state interactions (like, e.g., meson masses, 
decay constants, semileptonic form factors and e.m. transitions), the use of dif- 
ferent EC's on valence and sea quarks produce finite-volume effects which remain 
exponentially small. In this way there is no need for producing new gauge config- 
urations for each quark momentum, and this is quite relevant in the case of gauge 
configurations with dynamical fermions. 

The corresponding quark propagator S{x, y) = {iIj{x) ip{y)) still satisfies Eq. ([5]) 
with the same Dirac operator D{x,z) but with different EC's: 

^D(x,z) Siz,y)=5,,y (9) 

z 

Technically in order to work always with fields satisfying periodic EC's in space 
and time we follow Refs. [311 132] by introducing a new quark field as 

^l;~{x) = e-^^'^-^^'^ij^x) (10) 

where the four- vector 9 is given by {L/2T,6). In such a way the new quark 
propagator S^{x,y) = {'ipgix) V^gd/)) satisfies the equation 

J2d\x,z) S\z,y) = 6.,y (11) 

z 

with a modified Dirac operator D^{x, z) but periodic EC's in both space and time. 
The new Dirac operator is related to Eq. ((61) by a simple re-phasing of the gauge 
links 

U,{x) ^ Ul{x) = e2--^~-/^ U,ix) . (12) 

In terms of S^{x,y), related to the quark fields i'gix) with periodic EC's, the 
all-to-all quark propagator S{x,y), corresponding to the quark fields iplx) with 
twisted EC's, is simply given by 

S{x,y)=e''^''<^-yy'^ S\x,y) . (13) 

3 Two- and three-point correlation functions 

We are interested in the calculation of the vector form factor of a charged pion 
defined through the relation 

(7r+(pO|t>,(0)|vr+(p)) = F^{q') {p + p'), , (14) 
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where p {p') is the initial (final) pion four-momentum, = {p — p'Y is the squared 
four-momentum transfer and is a conserved e.m. current on the lattice. Split- 
ting into an isovector and an isoscalar part, it is easy to show that the matrix 
elements of the isoscalar component between pion states is vanishing in the con- 
tinuum limit, thanks to charge conjugation and isospin symmetric^. 

Thus, up to discretization effects we take at a generic (Euclidean) space-time 
point X = {tx, x) in the following form 

V,{x) = Zy V^{x) , (15) 
^m(^) = 2 [^(^)7M^(a;) - d{x)^^d{x)\ (16) 

with Zy being the renormalization constant of the isovector part of the vector 
current at maximal twist (cf. Ref. |28j). 

The insertion of the current ffTBl) generates two types of Feynmann diagrams, 
the so-called connected and disconnected diagram^l. In the former the external 
current is attached to the valence quarks, whereas in the latter the current inter- 
acts with the sea quarks. However, in the continuum limit the vanishing of the 
pion-to-pion matrix element of the isoscalar current {u'~f^u + d'y^d) implies that the 
connected diagrams stemming from the u- and d-quark terms of the current (fT6|) 
are equal in absolute value and opposite in sign, while their disconnected coun- 
terparts are vanishing for each quark flavor. Thus, in tmLQCD the disconnected 
diagrams for the e.m. pion form factor represent a pure discretization effect, which 
turns out to be of order O(a^) (see later on). Therefore, up to discretization ef- 
fects, it is enough to consider only the connected insertion of one single flavor of 
Eq. (dS]). 

From Eq. ffT^ the pion form factor can be extracted from both the time and the 
spatial components of the vector current. However for reasons which will become 
clear during this Section, we work in the Breit reference frame where p' = —p, so 
that the spatial components of the vector current are vanishing identically. There- 
fore we limit ourselves to consider the following two- and three-point correlation 
functions 

C-{t,p) = J2(0n{x) Oiiz)) 6,,..,^ e-^^^<^^-'\ (17) 

x,z 

Cr{t,t\p,p') = Y.^OMV,{x)Oi{z))5t,t^.t.5,,t,-t^ 

x,y,z 

^-ip-{x~z)+ip'-{x-y) Qg\ 



3t 



'Note that in tmLQCD the charge conjugation symmetry is preserved, while the isospin one 
is broken at finite lattice spacings. 

"^The terms "connected" and "disconnected" refer to fermionic lines only. 
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where Vo{x) = u{x)'you{x) and Ol{z) = u{z)'y^d{z) is the operator interpolating 
the TT"*" mesons. Note that, since we want to use all-to-all propagators, in Eqs. ( fT7|) 
and ( ITSl) there is an additional sum over the space-time lattice volume, which 
helps improving the signal quality with respect to the case of a fixed-point source 
(z = 0). 

Using the completeness relation and taking t and (t' — t) large enough, one gets 

(,-En{p)t ^-E^(p'){t'-t) ^ 1^20) 

where, up to discretization effects, Et,{p) = a/ + and y/Z^ = (0|O7r(0)|7r+) 
is independent on the meson momentum p. Note that in tmLQCD at maximal 
twist the value of the coupling constant determines the pion decay constant /^r 
[2H] without the need of the knowledge of any renormalization constant, namely 

A = , (21) 

TT 

where m is the bare twisted quark mass. 

Taking advantage of the choice of the Breit frame where p' = —p, it follows 

Crit,t',p,-p) 1 {n+ip')\Vo\n+{p)) 

where 

^ [EM - EAvt - \P- Pf - 4bf (23) 

is independent of the simulated pion mass. 

The vector renormalization constant can be obtained from Eq. (1221) by using the 
absolute normalization of the pion form factor at = 0, namely F^^^q^ = 0) = 1, 
which implies 

„-lT,-„ c;'(f,f',o,o) 

Combining Eqs. ( 122|) and ( l24l) one gets 

i?o(t, t'; q^) . 91M^ _ ^^(,.) (25) 

Q-(t,t',0,0) C-(t',p) (/_T)!:^ 
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which means that the pion form factor can be obtained directly from the plateau 
of the double ratio given by the l.h.s. of Eq. ( l25l) at large time distances. Note 
that in this way the normalization condition Fj^lq"^ = 0) = 1 is fulfilled at all quark 
masses, lattice volumes and spacings. 

The (mass-dependent) renormalization constant Zy can be obtained alterna- 
tively using the 3-point correlation function calculated in a frame in which the 
initial and final pions have the same momentum p, i.e. from the plateau of the 
ratio C'^{t',p) / C^'^it^t' .p^p) at large time distances. In this way the pion form 
factor can be extracted from the plateau of a ratio of 3-point correlation functions 
only, i.e. from C'^'^it.t' ,p,—p) / Cq'" (t , t' , p, p) . Such an alternative approach has 
been tested in Ref. and shown to have a statistical precision similar to the one 
based on Eq. ( l25l) . 

In terms of the all-to-all quark propagators S'^^'^\x,z), where the flavor labels 
u and d correspond to T3 = ±1 in Eq. (I6l), the two-point function C^it^p) of the 
charged pion becomes 

C^it.p) = J2{Tr[S.{x,zh,S,{z,xh,]) 6t,t^.t. e-'^^^^^-^. (26) 

x,z 

For the tmLQCD action the 75-hermiticity property 

Sd{z,x) =-f5Si{x,z)-f5 (27) 

holds with the dagger operator acting in the (suppressed) color and Dirac spaces. 

As for the three-point correlation function Cq'^ {tx^ty^p, —p), according to the 
discussion on the disconnected diagrams made before Eq. ( |T8l) , up to discretization 
effects one gets 

Co^'it^t' ,p,-p) = ^{Tr[Su{x,z)-f^J:duiz,x;t';-p)-fo]) 

x,z 

■ kt.-t. e-^'P<^-'^ , (28) 
where T,du{z, x; t'; -p) = 75[Sd„(x, z; t'; -p)]^75 and 

Erf„(x,^;t';p) = J2Sdix,yh,Suiy,z) e-''^<^^~y^ 5,,t,-t. . (29) 
y 

The segwentm/ propagator T,du{x, z;t';p) satisfies the equation 

J2M^^y) ^du{y,z;t'-p) = ^5Su{x,z) 5,,t._t. e'^-^^-^ . (30) 
y 

As it has been shown in Ref. [28], the calculation of correlation functions of 
parity symmetric operators is automatically 0{a) improved at maximal twist. 



11 



Thus for non-vanishing values of the spatial momenta the 0{a) terms can be 
eliminated by appropriate averaging of the correlation functions over initial and 
final momenta of opposite sign. Using the symmetry of correlation functions under 
the spatial inversion and the simultaneous exchange of u and d quarks (which is 
fulfilled only after gauge averaging at maximal twist in the physical basis) as well 
as the charge conjugation symmetry and the 75-hermiticity property, one gets 
that: i) the correlators (EE]) and ([28]) are real, and ii) C"'(V,p) = C'^{t\-p) 
and Cq'^IJ:, t',p, —p) = C^'^it, f, —p,p). Thus discretization effects in both C^it' .p) 
and Cq {t^t' ,p, —ff) start automatically at order (^(a^jfl. 

Let us now consider the case of quark fields with twisted BC's. Equations fl2Bl) - 
flHUl) hold as well by simply replacing the propagators S and S with the correspond- 
ing twisted ones, S and S, and by taking into account the change of the quantized 
momenta, namely pj —>■ pj [see Eq. (jl])]. The two- and three-point correlators 
can be expressed in terms of quark propagators satisfying periodic BC's, e.g. in 
terms of [see Eq. ( lT3i) ]. We now write down the explicit formulae for sake of 
completeness. 

In order to work in the Breit frame we consider three choices of the twisting 
four-vector 6, namely 6 = 6± = {L/2T,±6) and 6 = 6q = {L/2T,0) for various 
values of 6. Writing p in the generic form p = 27t6/L, we get 

C^it,-e) = J2{Tr[S:^{x,z)j,S'j>{z,x)^,]) (31) 

x,z 

Crit,t\^e,-^e) = J](Tr[^f+(a:,;^)754f"(^,^;07o])5M.-*. , (32) 

x,z 

where thanks to the 75-hermiticity property one has 

Sr-(^,x;0=75[S:-'^°]t(x,z;075 (33) 
and the sequential propagator S^~'^°(x, z; t') satisfies the modified Dirac equation 

5^D^(a:,y) T^J\y,z-t') = i.sl^x^z) . (34) 

y 

Note that, because of Eq. ffT^ . no exponential factors appear in the r.h.s. of 
Eqs. (13T1) -( 1M1) and the dependence on the vector 6 is totally embedded in the 
twisted quark propagators and S^-'^". 

^This result holds as well in all reference frames and it explains the findings shown in Fig. 12 of 
Ref. [16j , where the correlation functions with opposite momenta have been calculated explicitly 
and found to be identical within statistical errors. 
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3.1 Stochastic procedures 

The next point to be addressed is the evaluation of the all-to-all propagator S^{x, z) 
which is the solution of the modified Dirac equation (11 II) . Restoring color and spin 
indices, denoted by Latin and Greek letters respectively, one has 

Y}D\x, y)]% [S\y, z)]% = 5.,. 5,,, . (35) 
y 

The computation of exact all-to-all quark propagators is a formidable task well 
beyond present computational capabilities, because it involves a huge number of 
inversions of the Dirac equation for all possible locations of the source in space and 
time. Consequently most of the lattice computations of connected 2- and 3-point 
correlation functions are till now carried out using the point-to-all propagator by 
fixing the source at some space-time point, referred to as the origin. To get the 
expressions of our 2- and 3-point correlators in terms of point-to-all propagators 
it is enough to limit the sum over the variable 2; to ^ = everywhere in Eqs. ( l3Ti) - 
(IMl) . The basic advantage of the all-to-all propagator with respect to the point-to- 
all one relies in the fact that the former contains all the information on the gauge 
configuration, which in turn means that the calculation of 2- and 3-point functions 
using all-to-all propagators is expected to have much less gauge noise. 

An efficient way to estimate the all-to-all propagator is based on stochastic 
techniques with the help of variance reduction methods to better separate the 
signal from the noise (see Ref. ^7] and references therein). In recent years new 
stochastic methods have been developed, like the dilution method of Ref. [38j and 
the so-called "one-end-trick" of Ref. [30]. The latter, already applied by the ETM 
collaboration to the calculation of neutral meson masses (see Refs. [21] and [39j), 
allows to achieve a great reduction of the noise-to-signal ratio and it will be applied 
in this work to the calculation of 3-point correlation functions (see also Refs. [22] 
and [10]). 

The starting point of all stochastic approaches is to consider random sources 
77" (x), which, for reasons that will become clear later on, we take independent of 
both the spin variable and the twisting vector 9 (i.e., of the quark momentum). 
The index r (r = 1, ...N) enumerates the generated random sources, which must 
satisfy the following constraint 

1 ^ 

hm Y.r^^X^)[r,l{y)]* = 5a,, 5,^y . (36) 

N^oo iV •^"^ 
r=l 

In this work we adopt for the sources a random choice of ±1 values. Then one 
introduces the "0-propagator" 

y 
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which is solution of the equation 

Y,\D\x, y)]% [4iy%^ = ri^^x) S^,, . (38) 
y 

where the sum over repeated color or spin indices is understood. As explained 
in details in Ref. [21], the quantity X]r=i[*i^r(^)]a/3['7r(l/)]* is an unbiased 

estimator of the all-to-all propagator [S^{x,y)]'^^. However, while the signal is of 
order 0{1), the noise is of the order ^JV/N (where V is the space-time volume) 
and therefore a huge number of random sources and inversions of Eq. flH^ would 
be required. 

The "one-end-trick" is based on the observation that the product of two "0- 
propagators" is an unbiased estimator of the product of two all-to-all propagators 
summed over the intermediate space-time points. In this case, however, the signal 
is of order V, while the noise is of order V/\/N, so that it is even sufficient to 
employ one random source per gauge configuration, as we do in this work. 

Choosing the random source Tf{x) to be non- vanishing only for a randomly- 
chosen time slice, located at tr □, the 2-point correlation function (pTl) can be 
estimated as 

c'it^^e) = $^([0f^(x,t.)]^^ {K:MQr,ar kt.-u) m 

X,tx 

where we notice that the two 0's have the same flavor. Looking at the above 
equation the 0-propagator P^^Y^ ^ ^o\e quite similar to the one of the 

point-to-all propag ator [S%x, 0)]fi3 with only one color index, being the other one 
carried by the random source. This means that the time needed for the calculation 
of the 0-propagator is 1/3 of the one required for the point-to-all propagator. Note 
also that both 0^+(a;) and 0^°(a;) are solutions of Eq. (|38l) with the same random 
source rjr{x). This is essential to properly get the r.h.s. of Eq. flH^ . Moreover the 
independence of the random source from spin indices allows to evaluate 2-point 
correlation functions with interpolating flelds of the form {qVq') for any Dirac 
matrix F. 

The stochastic estimate of the 3-point correlation function ( l32l) requires the 
introduction of the sequential "^-propagator" 

y 

^The random choice of the time shoe at tr is mainly motivated by the reduction of autocor- 
relations observed for fermionic quantities using the ETM gauge ensembles (see Ref. ^Si)- 
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which is solution of the equation 



y 

One gets 

■ [757o]7« ^t,t.-u) ■ (42) 

Note that: i) the quark propagators required in Eqs. flHUl) and fH21) are those of 
one single flavor, while the other quark flavor appears only in the modified Dirac 
operator of Eq. ( HTl) . and ii) for each value of the quark momentum injected via 
the twisted BC's a new inversion of the Dirac operator is required. 



4 The charged pion form factor 

As already mentioned in the Introduction, the ETM collaboration has started an 
intensive, systematic program of calculations of three-point correlation functions 
relevant for the determination of meson form factors at low, intermediate and 
heavy quark masses. In this work we concentrate on the results obtained for the 
vector form factor of the pion. 

In Table[I]we collect the simulation set-up for all the runs carried out at /3 = 3.9 
and for the two runs performed at a finer lattice spacing {(3 = 4.05). Approximate 
values of the (charged) pion mass in physical units as well as of the quantity 
M^L, which governs finite volume effects in the so-called p-regime of ChPT, are 
reported for each run. 

The gauge configurations used for the measurements are selected from the 
trajectories produced by the ETM collaboration (see Refs. [23l |2ll |25]) at various 
values of the sea (bare) quark mass, arusea- We have chosen 1 configuration out of 
at least 20 (equilibrated) trajectories in all cases except for the run at the lightest 
pion mass (1 out of 10). 

Volume effects can be checked through the runs i?2a and i?26 (see later subsec- 
tion while lattice artifacts can be studied by means of the runs and R2C 
at a pion mass around 300 MeV and of the runs R^a and i?5ft for M^^ ~ 480 MeV 
(see later subsection 14. 5p . 

In this work at each value of the (bare) quark mass, am = arrisea, the statistical 
errors are evaluated with the jackknife procedure, while a bootstrap sampling will 
be applied in order to combine the jackknives for different quark masses (see later 
Section E]) . 
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R2b 
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Ba 


01 on 


24^ • 48 
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~ 5 1 


240 






Rq 


B^a,b 


0.0150 


24^ • 48 


~ 580 


~ 6.1 


240 


4.05 


~ 0.07 


R2c 


Cl 


0.0030 


32^ • 64 


~ 300 


~ 3.4 


240 






R5b 


C3 


0.0080 


32^ • 64 


~ 480 


~ 5.4 


240 



Table 1: Set-up of the lattice simulations for the various runs considered in this 
work. 



4.1 Vector renormalization constant Zy 

In tmLQCD tuned at maximal twist the constant Zy renormalizes both the isovec- 
tor part of the (local) e.m. current [see Eq. (fT6l) ] and the isovector off-diagonal 
components of the (local) axial current, e.g. A^{x) = J(x)7^75m(x). Therefore the 
renormalization constant Zy can be calculated in two ways. The first one is from 
Eq. fl24p , which makes use of 2- and 3-point correlation functions and is equivalent 
to fix the absolute normalization of the pion form factor, -F^(O) = 1. The second 
way is from the non-singlet axial Ward Identity {WI), which, up to discretization 
effects, in tmLQCD at maximal twist reads as [28j 

Zyd^A^ix) = 2amP,{x) (43) 

where am is the bare quark mass and P5(x) = c/(x)75m(x) is the bare pseudo-scalar 
density. The presence of bare operators in the r.h.s. of Eq. ( H3|) is due to the fact 
that at maximal twist the mass renormalization constant is equal to the inverse of 
the pseudo-scalar renormalization constant, i.e. Z^ = Zp^. At zero momentum it 
follows 

Zy = 2am (44) 

with A-{t,6) = j:.,am^) a(^)) kt.-t.- 

The results obtained for the ratio given by the r.h.s. of Eq. (|2^ . evaluated for 
all the runs ai jS = 3.9 and V ■ T = 24^ ■ 48 a^, are shown in Fig. [I](a). The time 
distance t' between the time slices of the source and the sink is fixed at t' = T/2, so 
that the 3-point correlation function ffT^ becomes antisymmetric with respect to 
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t = T /2 and it can be appropriately averaged to reduce the statistical fluctuations. 
Moreover for finite time extension T the 2-point correlation function C'^{t^p) is 
symmetric with respect to t = T/2, so that a second exponential e"'^'"*-^*-^"*^ 
appears in Eq. (fT9!) and a factor 1/2 has to be applied to the r.h.s of Eq. (!24|) . 

From the plateau region denoted by the vertical dotted lines in Fig. [T]^a) an 
estimate of the renormalization constant Zy can be obtained at each value of the 
bare quark mass. The results are reported in Fig. [11(b) and compared with the 
corresponding results obtained from the WI using Eq. ( 14^ (see Ref. [S]). Both 
methods exhibit an extremely high statistical precision of the order of 0.3 %. 
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Fig. 1: (a) Ratio of 2-point and 3-point correlation functions given by the r.h.s. of 
Eq. evaluated for t' = T/2 at /3 = 3.9 and V ■ T = 24^ ■ 48 a^, versus 

the (Euclidean) time t in lattice units, (b) The vector renormalization constant 
Zy as obtained at different values of the bare quark mass in lattice units. Open 
dots correspond to the values extracted from the plateau region denoted by the 
vertical dotted lines in (a). Open squares are the results obtained from the WI using 
Eq. ( [^^y R^f- I4^ - The solid and dashed lines are simple linear interpolations 
of the lattice points and the full markers denote the corresponding values at the 
chiral point. 

The quark mass dependence visible in Fig. [11(b) is a pure discretization effect 
and is different between the two methods. It appears to be linear in both cases, 
which is not in contradiction with the 0{a) improvement, since terms proportional 
to a^mAqcD niay be dominant with respect to terms proportional to a^rn^. 

The extrapolations to the chiral limit should therefore coincide, providing the 
value of the renormalization constant Zy., which is indeed defined in such a limit. 
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From Fig. [D^b) it can be seen that the values obtained by a simple linear fit 
at the chiral point coincide nicely within quite small statistical errors, namely 
Zv = 0.61088(14) from Eq. (El and Zy = 0.61076(19) from the WI [Eq. (HH)]. 



4.2 Momentum dependence of the 2-point correlation func- 
tion 



The 2-point correlation function ( l39l) has been calculated for various values of 
the twisting angle 6 chosen always in the symmetric form 6 = {6, 6, 9) with 
e = {0.0,0.11,0.19,0.27,0.35,0.44}. The time behavior of the effective mass (or 
logarithmic slope) aM^ffit), defined as 



aM,ff{t) = log 



C"(t,27rg/L) 
C^{t + a,2TTe/L) 



(45) 



is shown in Fig. [2] for two (representative) values of at /5 = 3.9 and V ■ T 
24^ ■ 48 
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Fig. 2: Effective mass of the pion ( [75| ) versus the (Euclidean) time distance in 
lattice units for ~ 300 MeV (a) and ^ 440 MeV (b) at p = 3.9 andV-T = 
24^-48 a'^. The twisting angle 9 is chosen in the symmetric form 9 = {9, 9, 9). The 
dots, squares, diamonds, triangles, the full dots and the full squares correspond to 
9 = {0.0, 0.11, 0.19, 0.27, 0.35, 0.44}, respectively. The dashed vertical line is drawn 
at t/a = 10, where the ground state starts to dominate. 



It can be seen that the statistical precision is remarkably high and it allows to 
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extract quite precisely the energy E.„{p) [see Eq. ( fTOll ] corresponding to the pion 
ground state, which starts to dominate from t/a = 10. 

The values obtained for the pion energy E.„^{p) are shown in Fig. [3] as a func- 
tion of the pion momentum given hjp = 2tt6/L, always at /9 = 3.9 and V ■ T = 
24^ -48 a^. The lattice points appear to be in remarkable agreement with the 
continuum-like dispersion relation Et,{p) = a/ M^{L) + \p\'^, where Mt,{L) is the 
charged pion mass at finite volume. We have also checked that, assuming the con- 
tinuum dispersion relation for the energy, all the 2-point correlation functions of 
moving pions can be simultaneously fitted very well together with the one at rest 
using a single momentum- independent matrix element [see Eq. (1191) ]. These 
findings clearly indicate that, at least for t/a > 10, where the ground state domi- 
nates, and for the pion momenta considered in this study, the discretization effects 
on the 2-point correlation functions C'^{t,2TT6/ L) are almost the same as those 
affecting the correlator at rest C"^(t, 0), which were investigated in Ref. [25] and 
found to be small. 
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Fig. 3: Squared pion energy E'^{p) in lattice units, obtained from the time plateaux 
of the effective mass shown in Fig. [H (hy choosing the time interval 10 < t/a < 21 ), 
versus the squared pion momentum p^ = 3(27r^^/L)^ in lattice units, for M.^ ^ 
300 MeV (a) and ~ 440 MeV (h) at (3 = 3.9 and V ■ T = 24^ ■ 48 The 
solid line is the continuum-like dispersion relation El{p) = M^{L) + while the 
dashed line in (a), which can be hardly distinguished from the solid one, represents 
the modified dispersion relation predicted by partially twisted and partially 
quenched ChPT at NLO elaborated in Ref. 



The use of twisted BC's is expected to produce finite volume corrections to 
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the continuum-like dispersion relation. Such corrections have been investigated 
in Ref. |12] using partially quenched ChPT at NLO. In the case of charged pion 
and adopting twisted BC's for one flavor only, the pion momentum p acquires an 
additive correction term K, namely 

El{p) = Ml{L) + {p + Kf (46) 

where the components of the vector K are given by {i ^ j ^ k) 

= -^7^ / dr — e~^(^)' e(r, 9,) 0(r, 9,) e(r, 9,) (47) 
V^Jtt^ Jo 

with Q{t,9) and Q{t,9) being the elliptic Jacobi function and its derivative. Ex- 
phcitly one has e(r, 9) = E^=-oo e-"("+^)' and e(r, 9) = E^=-oo(^ + ^)e-^("+^)'. 

We have evaluated Eq. (H71) for the run i?26, which has the smallest value of 
M^L (see Table [I]). The results are reported in Fig. [HJ^a) (dashed line) and they 
clearly indicate the smallness of the volume corrections to the pion momentum 
and therefore to the continuum dispersion relation expected at NLO. Thus finite 
size effects may be limited mainly to the pion mass and thus expected to be small 
(see Refs. [231 l2l] ). This is confirmed by the lattice results shown in Fig. HI where 
the pion energies obtained in case of the runs i?2a and R2b, which differs only for 
the lattice size, are compared. 



4.3 Momentum dependence of the pion form factor 

The advantage of calculating the pion form factor using all-to-all propagators, 
evaluated by the one-end-trick procedure with twisted BC's, with respect to the 
standard procedure based on point-to-all propagators with fixed sources and (spa- 
tially) periodic BC's is illustrated in Fig.O From the run i?2b we choose a different 
number of gauge configurations for the stochastic and non-stochastic procedures 
in order to get the same total computational tim^. 

Despite the more limited ensemble of gauge configurations the stochastic ap- 
proach provides a much better precision at the two lowest values of (a factor 
between ~ 2 and ~ 3). It also allows a very good determination of the form factor 
at the two highest values of considered in this study, where the procedure based 
on point-to-all propagators fails to give reliable signals even in the presence of a 
larger ensemble of gauge configurations. 

As discussed in the previous Section, the pion form factor FT^{q^) can be de- 
termined from the plateau of the ratio Ro(t,t'; g^), defined by Eq. fl^ . at large 

^Let us remind that the one-end-trick requires less computational time for a single inversion 
of the Dirac operator (a factor of about 1/3), but for each quark momentum a new inversion is 
needed by the use of twisted BC's. 
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Fig. 4: Squared pion 
energy El{p) in lattice 
units at ~ 300 MeV 
and (3 = 3.9 for the two 
runs i?2a md R2b, per- 
formed at the volumes V ■ 
T = . 48 (dots) 
andV - T = 32^ ■ 64 
(squares). The values of 
the twisting angle 9 are 
chosen in such a way that 
9/L has the same values 
in the two runs. The 
solid and dashed lines 
represent the continuum- 
like dispersion relation 
El{p) = Ml{L) + Ipf . 



time distances. The quality of the time plateaux is illustrated in Fig. [Ul while the 
momentum dependence of the extracted pion form factor F7r(^^) is shown in Fig. [7] 
for various values of M,r at (3 = 3.9 and V ■ T = 24^ -48 a^. We have checked 
that different choices of the time interval for the plateau region lead to values of 
F-Kio^) which are largely consistent within the statistical precision. The values of 
the pion form factor obtained for all the simulations of Table [T] are reported in the 
Appendix. 

In the whole range of values of both and the quark mass, considered in this 
work, our lattice data can be fitted very nicely using a simple pole ansatz 



1 



q^/M^ 



(4J 



pole 



as it is shown in Fig. [71 For comparison we also show the predictions of the Vector 
Meson Dominance (VMD) model, in which the parameter MpoZe is fixed at the 
value of the lightest vector-meson mass {Mymd) taken from Ref. [13]. The values 
obtained for Mp^ie by fitting our lattice points at /? = 3.9 are given in Table [21 

From Fig. [7| it can be seen that the VMD prediction, which considers the 
contribution of the lowest vector resonance only, is not exactly fulfilled, since Mpoie 
turns out to be systematically lower than Mymd- However such a comparison 
might be plagued by systematic uncertainties affecting the lattice determination 
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Fig. 5: Pion form factor F^^lq'^) versus in lattice units for a simulated pion 
mass of ~ 300 MeV. The full dots are the results obtained using twisted EC's 
in the Breit frame and the one-end-trick procedure for calculating the all-to-all 
propagators for an ensemble of 80 gauge configurations taken from the run i?2b- 
The open squares correspond to the results of the standard procedure based on point- 
to-all propagators with fixed sources for 120 gauge configurations of the run i?2fe- 
In this case spatially periodic EC's are applied in the frame where the final pion 
is at rest (p' = 0) and the momentum of the initial pion is given by p = 27r/L 
{(1, 0, 0), (1, 1, 0), (1, 1, 1), (2, 0, 0)}. At the two smallest values of q^ and for the 
ensemble of gauge configurations considered, only the stochastic procedure provides 
time plateaux of enough good quality to allow the extraction of the pion form factor. 



of the lightest vector-meson mass particularly at the lowest values of the pion mass 
(see Ref. 53)- Nevertheless, a simple extrapolation of Mpoie to the physical point, 
based on a polynomial fit in terms of quark masses (see later subsection I5.3p . yields 
0.713 ± 0.044 GeV, which is lower than the VMD prediction 



the value M^^f 



]{/rphys 



M„ = 0.776 GeV from PDG 



Note also that, defining the squared pole radius in terms of Eq. 



as 



pole 



pole 



dppole^q 

dq^ 



(49) 



=0 



the VMD model leads at the physical point to 



0.388 fm\ which 



underestimates by ~ 15% the (quite precise) experimental value of the squared 
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Fig. 6: Ratio Ro{t,t'; q^), defined by Eq. ( fl3]) . at t' = T/2 versus the time distance 
t in lattice units, for ^ 300 MeV (a), ~ 380 MeV (b), ~ 440 MeV 
(c), ~ 480 MeV (d) at 13 = 3.9 and V ■ T = 24^ ■ 48 The full 

dots, open squares, full diamonds, open diamonds and full squares correspond to 
a^q^ = —0.01, —0.03, —0.06, —0.10 and —0.16, respectively. The dashed vertical 
lines identify the region 10 < t/a < 14, where both the initial and the final pion 
ground states are isolated, so that the pion form factor F^.i^q'^) can be extracted. 



pion charge radius, {r'^Y^^' = 0.452 ± 0.011 fm^ [22]- On the contrary the value 
Ml^y^ = 0.713 ± 0.044 GeV imphes r^^;^ = 0.459 ± 0.057 fm^ in nice agreement 
with the experimental charge radius. 
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Fig. 7: Pion form factor F^(g^), extracted from the plateau region 10 < t/a < 14 
of the ratio RQ(t,T/2; q"^) (see Fig. IE), versus the squared 4.-momentum transfer 
q^ m lattice units, for ^ 300 MeV (a), =i 380 MeV (b), ~ 440 MeV 
(c), ~ 480 MeV (d) at /3 = 3.9 and V ■ T = 24^ -48 a^. The solid line is 
the pole behavior ( f^gp with the parameter Mp^ie fitted to the lattice points, while 
the dashed line is the VMD prediction with Mp^ie fixed at the value of the lightest 
vector-meson mass taken from Ref. l4^ . 



4.4 Finite size effects 

We have investigated the effects of the finite spatial extension L of our lattice boxes 
by comparing the results of runs i?2a and i?2b- In our simulations the latter has the 
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Run 


M (MeV) 


V -T / 


aM 1 


R^ 


~ 260 


32^ ■ 64 


359 ± 016 




~ 300 


32^ ■ 64 


363 ± on 


R2b 


~ 300 


24^ ■ 48 


379 ± 010 


R3 


~ 380 


243 ■ 48 


0.399 ±0.014 


i?4 


~ 440 


243 ■ 48 


0.420 ±0.012 


Rba 


~ 480 


243 ■ 48 


0.419 ±0.011 


Rq 


~ 580 


243 ■ 48 


0.440 ± 0.005 



Table 2: Values of the fit parameter Mp^e appearing in Eq. ([7^ obtained at j3 = 3.9 
in lattice units. 



smallest value of the quantity Mt^L, which governs finite size effects (FSE) in the 
p-regime. The physical extension of the two boxes is L ~ 2.8 fm and L ~ 2.1 fm, 
respectively. The values of the angle 6 are chosen differently at the two volumes 
in order to keep the values of fixed. 

The results for the pion form factor are shown in Fig. [HI while a direct compar- 
ison of the results for the pion mass and decay constant as well as for the squared 
pole radius, r^^i^, is illustrated in Table [31 
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Fig. 8: Pion form factor 
RiT^'f) obtained for the 
runs R2a (open squares) 
and R2b (full dots ), which 
correspond to different 
lattice boxes of size L ~ 
2.8 fm and L ^ 2.1 fm, 
respectively. The solid 
and dashed lines are the 
results of the pole fit 
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Run 


L ifm) 








R2a 


~ 2.8 


0.13377 (24) 


0.06625 (16) 


45.5 ±2.8 


R2b 


~ 2.1 


0.13623 (65) 


0.06459 (37) 


41.7 ±2.3 



Table 3: Values of the pion mass and decay constant from the high- statistics work 
of Ref. and of the squared pole radius [see Eq. (fTPP/ in lattice units for the 
runs i?2a o^nd R2b- 



It can clearly be seen that FSE effects are larger on the pion form factor (or, 
equivalently, on the pole radius) with respect to the case of the pion mass and 
decay constant. They indeed amount to ~ 8% on r^^^^ in contrast to a ~ 2% 
effect in the case of and /^r. However we notice that FSE effects on r^^^g are 
comparable to our statistical precision (~ 6%), while they are much larger in the 
case of and /^r (0.2-^0.6%). Thus it is mandatory to include volume corrections 
to our results at least on the pion mass and decay constant. 

On the theoretical side FSE on M^r and /j^ have been studied with ChPT at 
NLO in Ref. |44j and using a resummed asymptotic formula in Ref. |15], where 
both leading and subleading exponential terms are taken into account and the 
chiral expansion is applied to the vr — tt forward scattering amplitude. When the 
leading chiral representation of the latter is considered, the resummed approach 
coincides with the NLO result of Ref. [44] . Viceversa at NNLO the resummation 
technique includes only a part of the two-loop effects as well as of higher-loop 
effects. Recently the resummed approach has been positively checked against a 
full NNLO calculation of the pion mass in Ref. [IH], showing that the missing 
two-loop contributions are actually negligible. 

The volume corrections predicted by the resummed approach have been already 
considered in the analysis of the ETMC results for and /vr carried out in 
Refs. 1231 El [23 • 

On the contrary, till now, the theoretical investigation of FSE on the pion form 
factor is limited to the application of ChPT at NLO only. The case of periodic 
EC's is considered in Ref. [17], while twisted EC's are studied in Refs. [121 i3H] 
adopting two different reference frames, namely the rest frame of the final meson 
[i2j and the Ereit one [H] . 

The sign of the volume effects on the pion form factor depends crucially on the 
absolute value and the spatial direction of the twisting vector 9. The sign of FSE on 
the charge radius turns out to be opposite between the cases of periodic (Ref. [47]) 
and twisted (Refs. [12], HH]) EC's. When periodic EC's are used the extraction of 
the charge radius requires the use of the smallest available momentum, which is 
equal to 27r/L. Such a restriction is absent with twisted EC's and therefore volume 
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effects are different. 

Moreover the volume corrections depend on the reference frame: in the rest 
frame, besides the usual term related to the difference between the infinite volume 
loop integral and the sum over quantized momenta, there are two further contri- 
butions [12] arising from isospin and hypercubic invariance breakings generated by 
flavor-dependent twisted BC's. Such two terms are vanishing in the Breit frame 
as shown in Ref. [?8]. 

Only the results of Ref. [IE], in which both the twisted BC's and the Breit 
reference frame are considered, can be directly applied to our data. Thus one gets 



F,(g2;L)-F,(g2;oo) = ^ \ [ dx h,2 

Jtt Jo 



h/1 




(50) 



with g2 = -A{2'nd/Lf and 



h/2 



2710 



-.Ml 



27r3/2L2 



dr —= 

It 



.1=1 



(51) 



where B(r, 6') is defined after Eq. (|471) . The NLO volume corrections on the pion 
form factor expected for our run i?26 do not exceed half of the statistical error, and 
they are even smaller in the case of the run i?2a at the largest volume. The FSE's 
predicted by Eqs. fl50H51l) are quite small and have the same sign for all the choices 
of the twisting angle Q made in this work. The NLO corrections go to the right 
direction decreasing slightly the differences between the pion form factor obtained 
at the two box sizes. 

As for the squared pole radius, the shift with the lattice volume reported in 
Table [3] has the same sign expected from the volume correction fl50|) . However 
the FSE calculated at NLO for the run i?26 corresponds to an increase of ^ 3% 
only, that is almost a factor 3 less than the observed FSE (~ 8%). This suggests 
that higher-order chiral effects might be relevant on the pion form factor still for 
M^L ~ 3, although our statistical precision (~ 6%) does not exclude FSE's on 

•pole 



as small as the ones predicted at NLO by Eq. (!50|) . 



In the case of our runs -Ri and R2ai which correspond to iVf^L ~ 4, the NLO 
volume corrections on r^^^^ are expected to be ~ 1%. After multiplying such a 
value by a factor «i 3 in order to take into account conservatively higher-loop 
effects, the expected FSE remains well below the statistical precision. 

Therefore in this work we decide to analyze our form factor data using only 
simulations with M^^L > 4, which means in practice that the run R^h is excluded 
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from our analyses of the pion form factor. On the contrary in case of the pion mass 
and decay constant we keep the run in the set of fitted data, but the FSE's, 
calculated through the resummed asymptotic formula of Ref. |15] at the NNLO 
accuracy for the vr — tt forward scattering amplitude, will be taken into account 
(see Sections and E]) • 

4.5 Discretization effects 

We have investigated the impact of lattice artifacts on the pion form factor by 
considering the runs i?2c and R^h at the finer spacing a ~ 0.07 fm (see Table [T]). 
These runs correspond to pion masses equal to ~ 300 MeV and ~ 480 MeV, 
respectively, which are very similar to those of the runs i?26 and R^a at a ~ 0.09 
fm, while the physical lattice size is almost kept fixed {L ~ 2.1 fm). Our results 
are shown in Fig. [9] in terms of the Sommer parameter tq instead of the lattice 
spacing a. The ratio rg/a has been determined in the chiral limit at the two lattice 
spacings in Ref . [2S], obtaining rg/a = 5.22±0.02 at /9 = 3.9 andro/a = 6.61±0.03 
at (3 = 4.05. 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 



Fig. 9: Results of the pion form factor F.,;{q^) versus in units of the Sommer 
parameter tq, obtained for the runs i?2b (full dots) and i?2c (open squares) at M^- ~ 
300 MeV in (a), and for the runs R^a (full dots) and R^b (open squares) at M^^ ~ 
480 MeVzn (h). The physical lattice size is the same in both runs ~ 2.1 fm). 
The solid and dashed lines are the results of the pole fit ^B) , while the dotted line 
in (a) corrects the dashed one for the pion mass difference (see text). 

It can clearly be seen that the size of discretization effects is comparable to 
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the statistical error at both pion masses. At the lowest pion mass there is a slight 
mismatch between the values of M^ro corresponding to the runs R2b and i?2c- 
Using the ChPT formulae at NNLO evaluated in Ref. [H], which will be used in 
the next Sections, and adopting for the relevant LEC's the values given in Ref. [33] 
we have estimated the correction due to the pion mass difference and applied it 
to the results of the run i?2c (see dotted line in Fig. The correction is small, 
but reduces the impact of discretization effects, which now in terms of r'^^ie 
exceed ~ 5% at both pion masses. 

A more complete investigation of the scaling properties of the pion form factor, 
which requires the study of its mass dependence at two additional values of the 
lattice spacing, is needed and it is in progress. 

In the next Sections continuum ChPT will be applied to the chiral extrapolation 
of the results of the runs -Ri, i?2a, -R3, -R4, R^a and i?6, which correspond to a single 
lattice spacing (a ~ 0.09 fm) and to a pion mass range between ~ 260 MeV and 
~ 580 MeV with M^L > 4. The impact of lattice artifacts will be estimated by 
substituting the results of the runs and R^a with those of the runs R2C and 
i?56, respectively. 



5 Slope and curvature of the pion form factor 

The slope s and the curvature c of the pion form factor are defined from the 
expansion in 

F^{q^) = l + sq^ + cq^ + 0{q^) . (52) 
In terms of the pole ansatz (148|) the slope is given by 

1 

-*- pole (K'X\ 

Spole - j./r2 ~ c ' y^'^) 
^^^pole " 



while the curvature is constrained to be 



.2 \ 2 



We have therefore compared the slope and curvature obtained from the pole ansatz 
fHHl) with those of a simple cubic fit in q^ 

Fi'^'\q^) = 1 + s,ub + Cenfe + 4.6 ■ (55) 

The results obtained by including in the fitting procedure the form factor cor- 
responding to the four highest, negative values of are shown in Fig. [TU] (see also 
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Table Hj) □. It can be seen that the two determinations of the slope are in very good 
agreement and the results for the curvature are consistent within the statistical 
errors, which turn out to be lower in the case of the pole fit. 
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Fig. 10: The slope s 
(dots) and the curvature 
c (squares) of the pion 
form factor [see Eq. (5E) ] 
versus the squared pion 
mass in lattice units, for 
the runs Ri, i?2a; R3, Ri, 
R^a and Rq. Open dots 
and squares correspond 
to the results of the pole 
fit given by Eqs. ^53) and 
l^y, respectively. Full 
markers are the results 
obtained with the cubic fit 



(a M y 



In what follows we take as our best estimates the values of the slope and the 
curvature coming from the pole ansatz. The former ones, expressed in physical 
units using the value a = 0.087 fm from Ref. [23], are collected in the third column 
of Table H] and shown in Fig. [TT], where they are compared with the available re- 
sults of other lattice collaborations that employ 0(a)-improved lattice actions and 
unquenched gauge configurations. It can be seen that all the determinations of the 
pion charge radius exhibit a quite similar mass dependence, indicating that lattice 
artifacts are presumably under control. The results labeled as "QCDSF/UKQCD" 
in Fig. [TT]do not correspond to the original ones reported in Ref. p^. There the 
lattice spacing, instead of the Sommer parameter tq, was assumed to depend on 
the sea quark mass and such a procedure reintroduces non-negligible lattice arti- 
facts. The "QCDSF/UKQCD" results shown in Fig. [TT] are obtained after properly 

^The lattice points at the lowest, negative value of arc the noisiest data (see Fig.[7]and also 
Fig. [5] for the corresponding time plateaux) . The inclusion of these data in the fitting procedure 
does not change significantly the determination of the various parameters appearing in Eqs. (I53p 
and ((551). 
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extrapolating the ratio r^/a to the chiral limit§. 



Run 


(MeV) 


^polc — 6 s pole 

(fm^) 


^polc 

(10-3 fm^) 


1 — 6s^„K 

cub — ^"-'cub 

(fm^) 


^cub 

(10-3 fm^) 


Ri 


265 


0.352 ±0.030 


3.44 ±0.59 


0.364 ±0.028 


4.20 ± 1.12 


R2a 


304 


0.345 ±0.021 


3.30 ±0.40 


0.333 ±0.020 


2.45 ±0.82 


Rs 


383 


0.285 ±0.019 


2.25 ±0.31 


0.278 ±0.023 


1.85 ±0.54 


R4 


441 


0.258 ±0.015 


1.85 ±0.22 


0.249 ±0.021 


1.36 ±0.49 


Rba 


477 


0.259 ±0.014 


1.87 ±0.20 


0.254 ±0.018 


1.45 ± 0.44 


Rq 


584 


0.234 ±0.006 


1.53 ±0.07 


0.225 ±0.010 


1.14 ±0.26 



Table 4: Values of the pion mass, charge radius and curvature, determined from 
the pole ( f53l 54 ) and cubic ( f53]) fits, for the various ETMC runs. Physical units 
are used taking for the lattice spacing the value a = 0.087 fm from Ref. f23^ . The 
uncertainties are statistical (jackknife) errors. 



5.1 ChPT formulae at NNLO 



In Ref. [TT] the pion form factor, as well as the pion mass and decay constant, 
have been calculated in continuum SU(2) ChPT at NNLO in infinite volume using 
a modified minimal subtraction [MS) scheme to regulate the infinities. Using the 
quark mass rh as the expansion parameter, one has 



[^3 NNLO 



2Bm + K]nLO + [^']nNLO + ^(^' 



25m ■ 2x2 



2Bm -AxU — 
^ [N 

7 



2il + -L{f,) 



13 



± 



163 1 
96"iV2 



ki - 2k2 + 4£ 



3(^4 



9, 1, 

-k, + -k. 



A 



AI 



^-+2iV 



(56) 
(57) 



(58) 



fn 

[/^]nlo 



F + [AInlo + [AInnlo + ^( 



m 



2Fx2 



(59) 
(60) 



^We thank J. Zanotti for providing us the extrapolated values of rg/a at the chiral point. 
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Fig. 11: The squared pion charge radius versus the squared pion mass. Open dots: 
this work (see third column of Table Open squares: results from Ref. fWf 
corrected as explained in the text. Open diamonds, full triangles, full diamonds 



correspond to Refs. f21[ \2D^, respectively. The full dot represents the experi- 
mental value of the squared pion charge radius (r'^Y^'^' = 0.452 ± 0.011 frn^ from 
PDG 
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4.Fxl { ^ 
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-k4 + 3^6 
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(64) 



c = HnLO + HnNLO + ^(^) ' (65) 
2 

HnLO = .n/\rF4 > (66) 
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where 2Bm is the celebrated GMOR term, F is the pion decay constant in the 
chiral hmit (/^r is normahzed such that /^r ~ 130 MeV at the physical point) and 



N 


= (4vr)^ 






_ 25m 
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ki 






Am 


= 2£^ + 


2^(/^) ' 


A^ 




-L(^)] . 



(68) 

The constants are the finite part of the couphng constants appearing in 
the O(p^) Lagrangian after the apphcation of the MS procedure and their values 
depend on the renormalization scale /i through the anomalous dimensions 7^ as 
^'^M\/d^'^ = —'-)i/2N. The coefficients 7^ are calculated in Ref. [TU] and those 
relevant in this work are given by: 71 = 1/3, 72 = 2/3, 73 = —1/2, 74 = 2, 76 = 
— 1/3. The four constants r]^, r^, r[, r2 denote the contributions of the Oijf') 
Lagrangian after MS subtraction. Though the values of all the above constants 
depend on /i, at each order in the chiral expansion the physical observables are 
independent (as they should be) of the value of the renormalization scale /x. 
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At LO only two chiral parameters appear, namely B (related to the chiral 
condensate) and F. At NLO three further LEC's, £3, £4 and £g, are present. At 
NNLO the total number of LEC's increases up to 11 due to the inclusion of i^, i"^, 
^M, r}., rl and r^. 

We notice that the NNLO terms for the charge radius (1641) and the curvature 
(IU7I) do not depend upon the LEC's (i\ and separately, but only through the 
linear combination {i\ — ^2/2)- However different linear combinations of l\ and £2 
appear in the NNLO terms of both the pion mass (l58l) and decay constant (16T1) . 
Therefore the LEC's £^ and £2 can be determined by a simultaneous analysis of the 
the charge radius (and/or the curvature) together with the pion mass and decay 
constant. 

In what follows the 0{p'^) constants will be substituted by scale-invariant 
quantities, £i, defined via the relations 



The new quantities, which depend (logarithmically) on the quark mass, can be 
expressed as ii = log(A^/2i?m) and their values are commonly given at the physical 
point. 

We notice that in Ref. [TT] the quark mass is not actually used as the expansion 
parameter. Instead of it the physical pion mass and decay constant are adopted. 
In order to recover the formulae of Ref. [TT] it's enough to replace in Eqs. (I56|) - fl68|) 
X2 with M2//2, L(/i) with (1/A^) ■ \og{M^/fi'^) and to set Am = Ap = wherever 
they appear explicitly. 

As explained in Sec. 14.41 we apply to the pion mass and decay constant the 



corrections for FSE computed in Ref. |45j. Using again the quark mass m as the 
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expansion parameter, one gets 
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where -^^0,1,2 are modified Bessel functions, tlie values of the multiphcities m{n) 
are given in Ref. [15] and 
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Notice that in Eqs. (ffOl) - (ffTl) no further LEG is introduced with respect to Eqs. (l56l) - 

In order to recover the formulae of Ref. |45j , where the physical pion mass and 
decay constant are adopted as expansion parameters, it is enough to replace in 
Eqs. (I7n])-(I7ID X2 with M^//^, A^ with ^/n M^L and to set Am = Ap = 0. 



5.2 Chiral fits 

Let us now apply Eqs. (I56l)-( l67|) to the analyses of the quark mass dependence of 
our results. As explained in the previous Section the set of lattice data chosen for 
the fitting procedure is given by the results of the runs Ri, i?2a, R3, -R4, R^a and 
Re for four quantities: the pion mass and decay constant, the charge radius and 
the curvature of the pion form factor. In case of the pion mass and decay constant 
also the results of the run R2b are considered and the FSE corrections given by 
Eqs. (I70l)-(I7TD are applied. 

Since each run corresponds to an independent ensemble of gauge configurations 
a bootstrap procedure is applied in order to combine all the jackknives in different 
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ways (1000 samples are used in practice). The statistical uncertainties, which are 
reported here after, are therefore bootstrap errors. 

In order to fix the lattice spacing and the up/down quark mass the experi- 
mental values of the pion mass and decay constant (Mj?^^^- = 134.98 MeV and 
fphys. ^ ^3Q_7 ^ Q_4 j^gY f^Q^ 122]) are usecS- We determine firstly the 
value of the bare quark mass am^, at which the pion assumes its physical mass, 
by requiring that the ratio Mtt/Ztt from Eqs. (l56l) and (l59l) takes the experimental 
value 134.98/130.7 ~ 1.033. Secondly, using the physical value Z^^^**' the lattice 
spacing a is determined. The value of the renormalized light quark mass in the 
MS scheme, m'^^{2 GeV), is obtained from am.,^ by considering the determina- 
tion of the non-perturbative (multiplicative) renormalization constant = ^/Zp, 
evaluated using the RI-MOM scheme in Ref. [H], and the matching factor with 
the MS scheme, which is known up to four loops [5T] . 

We start with a ChPT analysis at NLO including our lattice data only for the 
pion mass and decay constant up to ~ 500 MeV in order to compare with 
the ETMC NLO analyses of Refs. [2Sl IMl ED- The main differences are: i) 
the use of a single lattice spacing (a ~ 0.09 fm) both in the present work and in 
Refs. [23l El], while the results obtained from two lattice spacings (a ~ 0.07 and 
0.09 fm) are taken into account in Ref. [3l]; ii) a better statistical accuracy of 
the data for M^, and /,r in Refs. [231 [211 El] due to the use of all the (correlated) 
trajectories produced by the ETM collaboration with respect to the present work 
in which only a subset of 240 (uncorrelated) trajectories are employed; iii) the 
presence of the results of the run Ri at M,r — 260 MeV in the present work and 
in Ref. [3l] at variance with Refs. p51[2l]. 

Following Ref. [21] the FSE correction can be evaluated beyond NLO using 
Eqs. fl 70 1) - 0711) and adopting for the unknown LEC's ii and £2 the central values 
given in Ref. [33j. The values obtained for the fitting parameters (the LEC's of the 
chiral Lagrangian) are given in the second column of Table [5], while the best fit at 
NLO, including the FSE corrections given by Eqs. (1701) -( |7T1) . is shown in Fig. [12] 
by the dashed lines. 

It can be seen that the values of all the chiral parameters, in particular of the 
LEC's £3 and £4, as well as the values of the lattice spacing a and the light quark 
mass mP'*2/s-^ ^re consistent with the findings of Refs. [23l [2^ Note also that 
the statistical precision of the extracted values of £3 and £4 are very similar in this 
work and in Refs. [231 [22 [Mj. 

The theoretical evaluation of FSE effects given by Eqs. (j7n|) -( !7T|) appears to 
work quite well. Indeed after applying the FSE corrections, the pion masses and 
decay constants corresponding to the runs i?2a (at L = 32a) and i?2fe (at L = 24a) 

^"^In order to account for the e.m. isospin breaking effects which are not introduced in the lattice 
simulations, we use the experimental value of the neutral pion mass in accord with Refs. [49|,l50j. 
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parameter 


NLO 


NNLO 


NNLO ± < r2 >^^P- 


2B (GeV) 


5.21 ±0.05 


5.19 ±0.42 


4.89 ±0.08 


F (MeV) 


121. 7± 1.1 


121 ± 10 


122.5 ±0.9 




3.48 ±0.12 


5.0 ±4.0 


3.1 ±0.7 




4.67 ±0.06 


5.2 ± 1.8 


4.39 ±0.14 


^6 


14.59 ±0.03 


14.9 ±2.2 


15.9 ± 1.3 




—0.4 [33j 


0.3 ± 3.4 


— 1.3 ± 1.4 




4.3 [33J 


5.3 ± 1.1 


5.1 ± 1.1 


M 




-0.1 ± 1.1 


-0.60 ± 0.29 


r'p ■ 10^ 




-0.5 ±1.5 


-0.15 ±0.14 


• 10^ 





-0.95 ±0.15 


-0.92 ±0.15 


rl ■ 10^ 





0.70 ±0.17 


0.77 ±0.04 


a (fm) 


0.0861 ±0.0007 


0.0863 ±0.0047 


0.0884 ±0.0006 




3.60 ±0.08 


3.6 ±0.6 


3.81 ±0.08 


/ d.o.f. 


0.72 


0.91 


0.88 



Table 5: Values of the chiral parameters, the lattice spacing and the renormalized 
quark mass ih = m^^^{2 GeV) at the physical point for various ChPT analyses 
(see text). For consistency with m, the parameter B is given in the MS scheme 
at a scale equal to 2 GeV. The values of the parameters r\.j, rp, r\ and r2 are 
given at the p-meson mass scale. In the case of the NLO analysis the parameters 
El and I2 O'fe. used only for evaluating FSE's, while the parameter Iq is fixed by 
the experimental value of the pion charge radius. The latter is not included in the 
NNLO analyses. The uncertainties are statistical (bootstrap) errors only. 



become consistent within one standard deviation, as shown in Table El 





lattice (%) 


theoretical (%) 


aM^{L = 2Aa)/aM^{L = 32a) - 1 


+ 1.8 (5) 


±1.2 


af^{L = 2Aa)/af^{L = 32a) - 1 


-2.5 (6) 


-2.5 



Table 6: Values of the quantities [aM.,^{L = 24a)/aM.,^{L = 32a) — 1] and [af^^lL = 
24a)/af-,r{L = 32a) — 1] from the runs md -R26 <^t a pion mass of c:^ 300 MeV. 
The lattice results correspond to the values given in Table\^ The theoretical results 
are those corresponding to Eqs. ([7^-([7l]j using for the relevant LEC's the values 
reported in the second column of Table\^ 

We have checked that the full exclusion of the run i?2b from our analyses does 
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Fig. 12: The ratio of the squared pion mass to the renormalized quark mass rh 



m^'^{2 GeV) (a) and the pion decay constant (h) versus m in physical units. The 
points at the highest value of m are not included in the analysis. The dots are 
the ETMC results, corrected by the FSE effects given by Eqs. ((7^-(|771j, and the 
squares represent the experimental value for each quantity from Ref. f2^ . The 
dashed lines correspond to the region selected at la level by the ChPT analysis at 
NLO [see Eqs. ( f5^) -([57y and ^5Mj-^63\)]. The values of the fitting parameters are 
listed in the second column of Table\^ 



not have any significant impact on the chiral fits as well as on the values obtained 
for the chiral parameters. Therefore, since the runs R2a and i?26 are compatible 
once theoretical FSE's are included through Eqs. (|70l)-( |7Tl) . in what follows we 
shall not show the results of the run i?26 in the figures (i.e. we show only lattice 
data having M^-L > 4) and we will always apply to the pion mass and decay 
constant the FSE corrections given by Eqs. ( !70l) -( !7Tl) . 

The quality of the NLO fit shown in Fig. [12] is quite remarkable, leaving ap- 
parently little room for higher-order corrections even at the highest pion mass 

580 MeV), though the latter point is not included in the fitting procedure. 
However we now show that the same does not hold for the charge radius and the 
curvature of the pion form factor. 

The NLO prediction for the charge radius [see Eq. (1631) ] depends in practice 
only on one LEG, Iq, being F fixed by the analysis of the pion decay constant with 
a precision of the level of ~ 1% (see the second column of Table [5]). Note also that 
both the derivative of [('"^)]nlo '^i^h respect to the quark mass and the curvature 
HnLO (^^^ ^'-l- independent of Iq and therefore basically parameter-free. 
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The value of the LEG ^Q can be determined from the experimental value of 
the squared charge radius, (r'^y^^' = 0.452 ± 0.011 fm^ [22], since the latter is 
expected to be dominated by the NLO term (163|) . This leads to 1% = 14.59 ± 

0. 03 (see the second column of Table [5]). The corresponding NLO predictions for 
the charge radius and the curvature are shown in Fig. [T3] by the dashed lines. 
They significantly overestimate our lattice data for the charge radius and largely 
underestimate those for the curvature. 

Alternatively we have excluded the experimental value of the charge radius 
and included in the fitting procedure the lattice data of the charge radius for pion 
masses up to ~ 500 MeV (i.e. for m < 0.05 GeV) obtaining Eq = 11.6 ± 0.3 (with 
X^/d.o.f. ~ 1.2). The corresponding NLO predictions are shown in Fig. [iniby the 
dotted lines. At the physical point the charge radius is {r'^Y^v^- = 0.352±0.008 fm 
in clear contradiction with the experimental valuj"1. 

Moreover, since the curvature is independent on i^, any NLO fit is unable 
to provide enough large values of the curvature consistent with the relation (13^ . 

1. e. with the pole behavior of the pion form factor observed in subsection 14.31 Such 
a finding suggests that a NLO analysis of the pion form factor is applicable only 
to quite low values of |g^| (see later Section [6]). 

Thus, both the NLO results shown in Fig. [13] and the smallness of the system- 
atic effects due to finite volumes and lattice spacings, estimated in subsections 14.41 
and 14. 5t indicate that the quark mass dependences of our lattice data for both the 
charge radius and the curvature require to take into account chiral effects beyond 
the NLO. 

The results of the fit performed using ChPT at NNLO [i.e. based on Eqs. (l56l) - 
(1671) ] are shown in Fig. [HI while the values of the fitting parameters are listed in 
the third column of Table [S] with the renormalization scale /z fixed at the physical 
p-meson mass. Notice that the experimental value of the pion charge radius is not 
included in the fitting procedure. 

As already found in Ref. [21], the inclusion of NNLO effects leads to quite 
large uncertainties in the values of all the LEC's, in particular both for the LEC's 
^3 and ^4 appearing at NLO and NNLO, and for the LEC's ii and £2 appearing 
only at NNLO. Nevertheless the uncertainties in the chiral fits shown in Fig. [TH 
are of the order of the statistical errors in the mass range of the lattice points. 
This means that the large uncertainties reported in the third column of Table [5] 
are strongly correlated; the effect of the variation of one fitting parameter can be 
always compensated by those generated by the variations of the other parameters. 
However when we extrapolate the chiral predictions for the pion mass and decay 

^^If only the lattice data of the charge radius for the two lowest pion masses {Mt^ < 300 MeV) 
are considered, the value of Iq becomes 12.8 ± 0.5 (with x^/d.o.f. ~ 0.65) and the predicted 
charge radius at the physical point is {r'^)P^y^- = 0.393 ± 0.017 fm^, which still deviates from the 
experimental value by three standard deviations. 
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Fig. 13: The squared charge radius (a) and the curvature (b) of the pion form 
factor versus the renormalized quark mass m in physical units. The dots are our 
lattice results and the square represents the experimental value of the squared charge 
radius f2^ . The dashed and dotted lines represent the region selected at la level 
by the ChPT predictions at NLO given by Eqs. ^.''^d / fg^) . In the case of the 
dashed lines the value of the LECIq is fixed by the experimental charge radius, while 
in the case of the dotted lines it is obtained by including in the fitting procedure 
our lattice data of the charge radius for pion masses up to ~ 500 MeV (i.e. for 
m < 0.05 GeV;. 



constant outside the mass range of the lattice data towards the chiral point we 
end up with rather large uncertainties as shown in Fig. [T31 

Such a situation is clearly unsatisfactory both for a precise extraction of the 
leg's and for the extrapolation to the physical point. The inclusion of the ex- 
perimental values of the pion charge radius and curvature in the set of fitted data 
can obviously reduce the uncertainties in the extraction of the chiral parameters, 
but in this way the predictive power of the chiral fits is lost. 

Thus we look for an observable which should be: i) unrelated to the vector 
form factor of the pion, ii) known experimentally and iii) whose chiral expansion 
at NLO contains one of the LEC's, let's say £3 or £4. In this way the experimental 
value of such an observable, expected to be dominated by the NLO contribution, 
can constrain sufficiently the range of the variability of one of the LEC's. In turn 
this could be beneficial to reduce the uncertainties of all our fitting parameters. 

A possible, appropriate choice is the squared radius {r'^)s of the pion "scalar" 
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Fig. 14: The ratio M^/rh (a), the pion decay constant (b), the charge radius (c) 
and the curvature (d) of the pion form factor versus the renormalized quark mass 
171 in physical units. The dots are our lattice results and the squares represent the 
corresponding experimental values from PDG W^ . The dashed lines correspond to 
the region selected at la level by the ChPT fit at NNLO based on Eqs. ( f5^) -([g?l). 
The values of the fitting parameters are listed in the third column of Table S The 
experimental value of the pion charge radius is not included in the fitting procedure. 



form factor, defined as 



'dFlil 



(73) 
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where 



(74) 



Indeed, on one hand side the experimental value of the pion scalar radius is known 
quite accurately from the analysis of vr — vr scattering data (see Ref. [53]), which 
gives 



(r2)^-P- = 0.61 ± 0.04 fm' 



(75) 



On the other hand side the chiral expansion of {r'^)s, calculated at NNLO in 
Ref. [11], reads as 



ir')s = [(O.InLO NNLO + ^(^') 



(76) 



2 

[(^')5]nlo = 




-6L(.)-i|) , 



(77) 
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N 



AI 



(7J 



It can be seen that the LEG £4, which also governs the NLO correction to the pion 
decay constant [see Eq. (!60l) ]. appears in Eq. (1771) . 

As already stressed, the experimental value {r'^Yg^' is expected to be dominated 
by the NLO contribution (1771) . Using the values of the relevant LEC's of the second 
column of Table [5] one gets [{r'^)s\NLO = 0.716 ± 0.014 fm^, which overestimates 
the experimental value (1731) by almost three standard deviations. Thus we want to 
use the NNLO calculation of (r^)^, also for consistency with the use of Eqs. fl3Ul) - 
( 1671) for the other observables. In order to do that we need to set the value of the 
parameter appearing in Eq. (1781) . In Ref. [H] an estimate of at the p- meson 
mass scale has been obtained using a resonance model, namely —0.3-10"^. We 
have checked that using the above value or putting the parameter equal to zero 
does not produce any significant difference in our chiral fits. This is not surprising 
since the effects of a non- vanishing value of are expected to be relevant at large 
quark masses only. Thus in what follows the value = is assumed. 

Including the experimental value {t'^Y^^' in the ensemble of fitted data and 
Eqs. (I7U|) - (I75|) in the fitting procedure, we obtain the results shown in Fig. [T^ with 
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Fig. 15: As in Fig. [17} but including the experimental value of the pion scalar 
radius from Ref. fS^ in the fitting procedure. The resulting values of the fitting 
parameters are listed in the fourth column of Table\B 



the values of the fitting parameters reported in the fourth column of Table [5l Let 
us remind that the experimental value of the pion charge radius is not included in 
the fitting procedure. 

Our expectation about the reduction of the uncertainties of the fitting param- 
eters is fully confirmed. Thanks to the introduction of the experimental value 
^^2yxp. value of £4 is determined quite accurately and this is beneficial for re- 
ducing the uncertainties of all the other LEC's (compare the third and the fourth 
columns of Table [5]) . 
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Note that with respect to the NLO analysis the values of both the parameter 
B and the lattice spacing, obtained in the NNLO analysis which includes (t"^) tT^', 
change beyond the corresponding statistical errors. On the contrary the values of 
the LEC's F, £3, £4 and Iq do not change significantly. 

Since our data used for the curvature rely on the assumption of the monopole 
behavior (see Eq. fl54p ) . we have checked that the values of the LEC's, extracted by 
including the curvature data obtained from the cubic fit ( |55l) . change only slightly 
within the statistical errors with respect to the ones reported in the fourth column 
of Table El 

The values of the charge radius and the curvature predicted at the physical 
point by the chiral fit shown in Fig. [T51 are (r 

^phys. _ (^^11 -|- 0.47) ■ 10"'^ fm . We perform a rough estimate of the systematic 
errors due to finite volume and discretization effects. Firstly we substitute the 
run i?2a with the run i?2fe and fit the new set of data; the changes in the central 
values of the chiral parameters provide an estimate of the finite volume effects. 
Secondly we further substitute the runs and R^a with the runs R2C and R^b 
at the finer lattice spacing, respectively, obtaining an estimate of discretization 
effects. Adding all the systematic uncertainties in quadrature, our final results are 

(^^2yhys. ^ g_45g ^ g ggg ^ g g24 ^ ^jq-^ 

^Phys. ^ (5.11 ± 0.47 ± 0.41) ■ 10~3 fm^ . (80) 

where the first error is statistical and the second one systematic. We remind that 
the findings (I7U|) and flHUl) are not completely independent, because they are based 
on an ensemble of fitted data which satisfy Eq. ( 15^ . 

The results for the pion charge radius at the physical point obtained by various 
lattice collaborations performing unquenched calculations are compared in Table 
Hand in Fig. [M 



collaboration 


Nf 


Action 


Y -T 1 a'^ 


a (fm) 


AU (McV) 


(r'')P'^y=- (fm^) 


ETM [this work] 


2 


tlSym + tmW 


32^ ■ 64 


~ 0.09 


> 260 


0.456 ± 0.038 


JLQCD 20 


2 


Iw. + overlap 


16^ ■ 32 


~ 0.12 


> 290 


0.404 ± 0.031 


JLQCD 18 


2 


plaq. + Clover 


20^ ■ 48 


~ 0.09 


> 550 


0.396 + 0.010 


QCDSF/UKQCD 19 


2 


plaq. + Clover 


24^ ■ 48 


~ 0.08 


> 400 


0.441 + 0.019 


UKQCD/RBC [21] 


2 + 1 


DWF 


16^ ■ 32 


~ 0.12 


> 330 


0.418 + 0.031 


LHP 14] 


2 + 1 


Asqtad + DWF 


20^ ■ 64 


~ 0.12 


> 320 


0.310 + 0.046 



Table 7: Results for the pion charge radius extrapolated at the physical point by 
various lattice collaborations performing unquenched calculations. The uncertainty 
of the ETM result corresponds to the statistical and the systematic errors, given 
by Eq. [7^ , added in quadrature. 
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Fig. 16: Results for the 
pion charge radius ex- 
trapolated at the physi- 
cal point by various lat- 
tice collaborations per- 
forming unquenched cal- 
culations (see Table [^. 
The vertical dashed lines 
show the experimental 
value {r^y^P- = 0.452 ± 
0.011 frn^ 



Our finding ( 1791) agrees very well with the experimental value (^r'^y^P- = 0.452 ± 
0.011 fm^ [22]. It is also consistent within the errors with the results of JLQCD 
{Nf = 2), QCDSF/UKQCD {Nf = 2) and UKQCD/RBC (A^^ = 2 + 1) collabo- 
rations, while the difference with the result of LHP {Nf = 2 + 1) collaboration is 
equal to ~ 2 standard deviations. 

We stress that the use of twisted BC's and the inclusion of the NNLO terms in 
the ChPT analyses are two important features considered in this work. Note that 
twisted BC's are used only in Ref. [21] and a NNLO ChPT analysis is carried out 
only in Ref. [20]. 

The ETMC result ( !79l) has been obtained by fitting the lattice data for values 
of the (squared) four-momentum transfer = —q^ up to 0.5 GeV^ using the pole 
ansatz ( 1481) or equivalently the monopole functional form 

p{^onapole),Q2^ ^ L . (81) 

1 + 

An interesting question is at which value of the predictions based on Eq. (ISTj) 
start to deviate from the experimental data. Recently, thanks to the CEBAF 
facility at JLab, the pion form factor has been measured quite accurately up to few 
GeV^ The comparison with the monopole prediction (IHTj) . using for the (squared) 
pion charge radius either the ETMC result (I79p or its experimental value from 
PDG [22], is illustrated in Fig. HH 

It can clearly be seen that there is no hint of a deviation of the experimental 
data from the monopole ansatz up to ~ 2 -h 3 GeV^. 

New experimental data up to — 6 GeV^, expected to be taken after the 
completion of the JLab upgrade to 12 GeV [S2], may shed light on the range of 
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Fig. 17: Pion form fac- 
tor times = —q^, 
Q'^Ft,{Q'^), versus in 
physical units. The 
dots, squares and dia- 
monds are experimental 
data from Refs. |3 0/ 
and 13 [3^, respec- 
tively. The dashed and 
solid lines correspond to 
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level by the predictions of 
the monopole form l[8l\) 
using (r^) = (j-'^y^^' = 
0.452±0.011 fm^ fMBJ and 
(r2) = {r^)phys. = o.456± 

0.038 fm^, respectively. 



validity of the monopole ansatz. 



5.3 Polynomial fit 

We want to discuss briefly an alternative flt to our lattice data based on a simple 
polynomial form, which does not have any logarithmic term; namely 



Ml 

TT 



2Bm ■ [1 + Cim + Cam^] 
F ■[! + Dim + D2rh^] , 
6/(Mo + Eirh + Eam^)^ . 



(82) 
(83) 
(84) 



The results of such a flt, applied to our lattice data having M^^L > 4 without 
applying any FSE correction, are illustrated in Figs. [18] and [T9l 

Among the fltting parameters we obtain 2B = 4.79 ± 0.08 GeV and F = 
126 ±2 MeV, while the lattice spacing turns out to be a = 0.0889 ± 0.0010 fm and 
the (renormalized) up/down quark mass is mP^^^' = 3.8 ±0.1 MeV. 

It can clearly be seen that the quality of the polynomial flt is quite similar to 
the one of the ChPT flt shown in Fig. [15] and therefore we have to conclude that 
our lattice results do not show a clear-cut evidence of chiral logs. 
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Fig. 18: The ratio M^/rh (a) and the pion decay constant (h) versus renormalized 
quark mass rh in physical units. The dots are the ETMC results, uncorrected for 
FSE effects, and the squares represent the experimental value for each quantity 
from PDG f2^ . The dotted lines correspond to the region selected at la level by 
the polynomial fits of the squared pion mass and decay constant ^8^. 

6 Chiral fits of the pion form factor 

In this Section we present the ChPT analysis of our results for the pion form 
factor including its momentum dependence. The ChPT expansion of FT^{q^) has 
been calculated in Ref. [TT] at NNLO, obtaining 



^-(^') = l+[^-(^')]l-loop+[^'^(^')]2-loops + ^(^' 



1-loop 



2-loops 



2xo 



-(^-4)JH-^(£^ + -L(/i) + - 



^5) 
^6) 



Axl Y'viw) + Uv{w) - AmJ{w) 



(87) 
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Fig. 19: The charge ra- 
dius of the pion form fac- 
tor versus the renormal- 
ized quark mass m in 
physical units. The dots 
are the ETMC results 
and the squares repre- 
sent the experimental va- 
lue of the pion charge ra- 
dius J2^ . which is not in- 
cluded in the fitting pro- 
cedure. The dotted lines 
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where w = /2Brh. The polynomial part Pv{w) is given by 
Pviw) 



w 



—ki H — /eo kA H — ka + r\ 

24 12 2^ 
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while the dispersive part Uyiw) reads as 
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where 



J{w) 


= z h{z) + — , 


(ill') 


— 7 h'^(7\ 


H2{w) 


iV 




= N-h\z) + —- 
w Nw 




= —(lH,iw) + 
wz \2 


1-A/w 


and 



(90) 
(91) 
(92) 



M^) - ^ , (93) 
ii^3(«^) + lj» + ^.) , (94) 

with z 

Using the above formulae it is possible to test the momentum dependence of 
the pion form factor predicted by ChPT at NNLO. Such a dependence is analytical 
up to the inelastic threshold g^^^ = 4M^. Thus, in the chiral limit, terms of the 
form log(— g^) appear in the pion form factor, which becomes a non- analytic 
function of q^. This is the origin of the divergency of both the charge radius and 
the curvature in the chiral limit [see Eqs. flU21) - flU7j) ]. 

It is easy to check that an expansion of Eqs. (!85l)-( !87|) in powers of leads 
to the result: FT,{q^) = 1 + (r^) g^/6 + c g^ + C(g^), with (r^) and c given by 
Eqs. (I62D-(I61 and Eqs. (I65l)-(I67D, respectively. Thus by using Eqs. (IH5l)-(IHZD it is 
possible to take into account (at least partially) the effects of order O(g^) in the 
momentum dependence of the pion form factor. 

Note that the NNLO terms fl871l89l) do not depend upon the LEC's and 
separately, but only through the linear combination {l\ — ^2/2)- Since different 
linear combinations of i\ and £2 appear in the NNLO terms of both the pion 
mass fl58|) and decay constant fl6T]) . the LEC's (i\ and £3 can be determined by a 
simultaneous analysis of the form factor together with the pion mass and decay 
constant. 

As already discussed in subsection 15.21 at NLO the pion form factor depends 
on one LEG, -^e, which governs only the linear term in g^. The ChPT predictions 
at NLO corresponding respectively to 1^ = 14.59 ± 0.03 and Iq = 11.6 ± 0.3 (with 
2B = 5.21 ± 0.05 GeV and F = 121.7 ± 1.1 MeV), reported already in Fig. [TS]in 
the case of the pion charge radius and curvature, are shown in Fig. [20] for various 
values of the pion mass. We remind that the value = 14.59 ± 0.03 is fixed by 
the reproduction of the experimental charge radius, while the value = 11.6 ±0.3 
is obtained by fitting our lattice data of the charge radius for pion masses up to 
~ 500 MeV. It can be seen that: 
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• the momentum dependence predicted by ChPT at NLO is almost linear at 
variance with the pole behavior f HHl) observed in our lattice data [see also 

Fig.Mh)] ; 

• using = 11.6 ± 0.3 the NLO approximation appears to work up to = 
—q^ ~ 0.15 GeV^ and for pion masses below ~ 300 MeV. With such a value 
of 4 the NLO formula ([63]) yields {r^)P^'y'- = 0.352 ±0.008 fm^ which under- 
estimates significantly the experimental charge radius (see subsection 15. 2p . 
A slight improvement can be achieved by using directly the NLO formula 
f lHUj) to fit our lattice data for the pion form factor at the lowest Q^-value (~ 
0.05 GeV^) and for the two lowest pion masses (M^ ~ 260 and ~ 300 MeV). 
We obtain 4 = 12.2 ± 0.5 corresponding to {r^)P''y'- = 0.373 ± 0.017 fm^ 
which still deviates from the experimental value by four standard deviations; 

• using 4 = 14.59 ± 0.03, which instead reproduces the experimental value of 
the pion charge radius, the range of applicability of the NLO approximation 
reduces to values of at least not larger than ^ 0.03 GeV^ and to pion 
masses below ^ 300 MeV, which are not covered by our present lattice 
data {Q^ > 0.05 GeV^). We notice that within the above restricted range of 
values of and pion masses the deviation of the pion form factor from unity 
becomes smaller than few percent and therefore a particular attention should 
be paid to the statistical precision as well as to the systematic uncertainties 
related to cut-off and finite size effects. 

From Fig. [20] it is clear that the description of our lattice data requires the 
inclusion of higher-order ChPT effects, which should be much larger than the 
finite volume corrections and the scaling violations observed in Figs. [H] and M for 
the values of and M^^ considered in our calculations. 

We have therefore performed a simultaneous NNLO fit of the lattice data of 
the runs Ri, i?2a, R3, -R4, -Rsa and -Re for the quantities M^, and Ft,{Q^), 
including all values of = —q^ from ^ 0.05 GeV^ up to ~ 0.8 GeV^. As in the 
previous Section, the constraint f[731) on the pion scalar radius is included in the 
fitting procedure in order to reduce the uncertainties in the extracted values of the 
chiral parameters. The latter are given in the second column of Table [H] The nice 
quality of the NNLO fit is illustrated in Figs. [2T]and[22] The corresponding value 
of the pion charge radius, calculated at the physical point using Eqs. (I62ti64p . is 
(^j.'iyphys. _ Q -|- 0.029 fm^, in nice agreement with the finding fl79|) . which, we 
remind, is based on the use of the pole ansatz fl48p that describes very well the 
momentum dependence of our lattice data (see Figs. 171 and 1201) . 

The comparison of the results shown in the fourth column of Table [5] and 
those in the second column of Table [8] clearly indicates that within the statistical 
uncertainties the extracted values of the chiral parameters are quite stable against 
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Fig. 20: Pion form factor Fj^{Q'^) versus the squared ^-i^omentum transfer = 
—q^ in physical units, for the run Ri at M^^ ~ 260 MeV (a), the run R2a o>t M^- ~ 
300 MeV (h), the run R3 at ~ 380 MeV (c) and the run R^ at ^ 440 MeV 
(d). For the lattice spacing the value a = 0.0861 fm is adopted from the second 
column of Tahle\^ The solid line is the pole behavior ( f^gp with the parameter Mpde 
fitted to the lattice points. The dashed and dotted lines are the regions selected 
at la level by the ChPT predictions at NLO corresponding respectively to = 
14.59±0.03 andie = 11.6±0.3 with2B = 5.21±0.05 GeV and F = 121.7±1.1 MeV 
(see text). 



chiral effects of order 0{q^) in the pion form factor with the only exception of the 
parameter which instead exhibit a rather large variation. The latter will be 
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parameter 


g2 < 0.8 GeV" 


g2 < 0.5 GeV" 


g2 < 0.3 GeV" 


2B (GeV) 


4.89 ±0.10 


4.91 ±0.07 


4.90 ±0.09 


F (MeV) 


122.6 ± 1.1 


122.5 ±0.8 


122.5 ± 1.0 




3.14 ± 1.03 


3.24 ±0.53 


3.19 ±0.81 




4.37 ±0.27 


4.41 ±0.10 


4.39 ±0.19 




15.0 ±0.9 


14.8 ± 1.1 


14.8 ± 1.5 




—0 54 ±101 


—0 31 ± 1 12 


— 28 ± 1 91 




4 40 ± 78 


4 29 ± 1 00 


4 23 ± 1 52 


r\, ■ 10^ 


-0.48 ± 0.30 


-0.43 ± 0.23 


-0.44 ±0.37 


■ 10^ 


0.11 ±0.19 


0.08 ±0.10 


0.08 ±0.20 


r\ ■ W 


-0.98 ±0.12 


-0.94 ±0.13 


-0.90 ±0.14 


rl ■ W 


0.43 ±0.03 


0.46 ±0.03 


0.52 ±0.05 


a (fm) 


0.0883 ±0.0006 


0.0883 ± 0.0006 


0.0883 ±0.0007 




3.80 ±0.09 


3.79 ±0.07 


3.79 ±0.09 


/ d.o.f. 


29 / 34 


19 / 28 


13 / 22 



Table 8: Values of the chiral parameters, the lattice spacing and the renormalized 
light quark mass m = m^^'^(2 GeV) at the physical point, obtained from the si- 
multaneous ChPT analysis of the pion mass, decay constant and form factor made 
at NNLO using Eqs. ( f5^) -( fg7]) and ( fg3]) -([g7|j, including also the constraint [75^ 
on the pion scalar radius. The second, third and fourth columns correspond to 
different Q'^-ranges of the lattice data of the form factor considered in the fitting 
procedure. The values of the parameters r\,^, rp, r[ andr^ are given at the p-meson 
mass scale. The uncertainties are statistical (bootstrap) errors only. 



included in the systematic error (see later Table [9]), providing the dominant source 
of uncertainty for the parameter rl^. 

We remind that, at variance with the results reported in Table El the ones 
shown in Table [HI are obtained without the assumption of the pole ansatz (HSl) for 
the momentum dependence of the form factor, but using only the functional forms 
(IHMHnD predicted by ChPT at NNLO. 

Effects from higher orders in the chiral expansion are expected to become more 
and more important as the value of increases. In order to check their relevance 
in our analysis we repeat the NNLO fit by limiting the range of values of Q^, i.e. by 
including only lattice data with < 0.5 GeV^ (see third column of Table [8]) and 

< 0.3 GeV^ (see fourth column of Table [H]). It can clearly be seen that, within 
the statistical precision, the extracted values of all the chiral parameters are only 
slightly sensitive to the Q^-range used and therefore to higher-order effects. 
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Fig. 21: The ratio of the squared pion mass to the renormalized quark mass rh 
(a) and the pion decay constant (b) versus m in physical units. The dots are the 
ETMC results and the squares represent the experimental value for each quantity 
from Ref. W^ . The dashed lines correspond to the region selected at la level by 
the NNLO ChPT analysis of the ETMC results for the pion mass, decay constant 
and e.m. form factor. The experimental value of the pion scalar radius ^J^j is 
added to the fitting procedure employing Eqs. ((76|j-([ygp. The values of the fitting 
parameters are listed in the second column of Table\^ 



Before closing this Section we mention that in Refs. [211 EHl [3lj a different 
definition of the LEC's at NNLO is adopted, namely the constants r]^ and r^. are 
replaced by the constants /cm and kp- Using the results of the second column of 
Table [H] we obtain fc^ = —0.6 ± 1.6 and kp = 1.2 ± 1.5. 



7 Final results for the LEC's 

In this Section we provide the final estimates of the LEC's from the present work. 
Our results, including both the statistical and the systematic uncertainties, are 
collected in the second column of Table O They have been evaluated by averaging 
the three central values reported in Table 8, using the quoted errors as the weights. 

As in the case of the charge radius and curvature, discussed in the previous Sec- 
tion, we estimate the systematic errors due to both finite volume and discretization 
effects. Firstly we substitute the run i?2a with the run i?2b and fit the new set of 
data; the changes in the central values of the chiral parameters provide an estimate 
of the finite volume effects. Secondly we further substitute the runs R2b and R^a 
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a 



1.0 



0.9 - 



0.8 



- " ^ ^ 



0.7 - 



0.6 
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0.00 



= 0.05 GeV^ 
= 0.15 GeV^ 
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= = = = = =C 



::;:0 




with <r^> 

s 

J I I I I I I I I I I L. 



0.02 0.04 0.06 
(GeV) 



0.08 



A 

m 



Fig. 22: ETMC results 
for the pion e.m. form 
factor versus the renor- 
malized quark mass m at 
various values of Q"^ (= 
—q^). The various lines 
correspond to the regions 
selected at la level by 
the ChPT fit at NNLO 
based on Eqs. 



with the fitting parame- 
ters given in the second 
column of Table The 
experimental value of the 
pion scalar radius [75\ ) is 
added to the fitting pro- 
cedure using Eqs. /fTT 



with the runs R2c and R^b at the finer lattice spacing, respectively, obtaining an 
estimate of discretization effects. All the systematic errors, which include also the 
spread of the central values of Table [HI are finally added in quadrature. 

In Table our estimates of the chiral parameters are compared with available 
results from NNLO ChPT analyses of vr — vr scattering data from Ref. [33], and 
with estimates obtained using VMD models in Ref. [H]. Our values for the LEC's 
^i, ^2, ^3 and £4 agree nicely with those extracted in Ref. [33]. The uncertainties 
obtained in this work for the LEG £4 is quite similar to the one from Ref. 
while ii and £2 are determined more precisely in Ref. [23] and £3 in the present 
work. 

On the contrary the estimates of the counter-terms r[ and r2 obtained in 
Ref. [11] adopting VMD models turn out to be much larger than our values by a 
factor ^2^3. 
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LEC 


ETMC (NNLO) 
(this work) 


non-lattice 
estimates 


2B (GeV) 


4.90 ±0.09 ±0.20 




F (MeV) 


122.5 ± 1.0 ± 1.0 






—0 4 ± 1 3 ± 6 


—0 4 ± 6 331 




4.3 ± 1.1 ± 0.4 


4.3 ± 0.1 [33] 




3.2 ±0.8 ±0.2 


2.9 ± 2.4 [33J 




4.4 ±0.2 ±0.1 


4.4 ± 0.2 1331 


^6 


14.9 ± 1.2 ±0.7 






-0.45 ±0.30 ±0.10 




■ 10^ 


0.08 ±0.16 ±0.05 




rl ■ 10^ 


-0.94 ±0.13 ±0.10 


-2.0 dH 


■ 10^ 


0.46 ±0.03 ±0.31 


2.1 [llj 



Table 9: Values of the LEC's obtained from the NNLO ChPT analyses of the 
previous Section and compared with available estimates arising either from NNLO 
ChPT analyses o/vr — vr scattering data fS^ or from VMD models 111]. The values 
of the parameters r\j, rp, rl and r2 are given at the p-meson mass scale. In the 
second column the first error is statistical and the second one systematic. 



The results obtained for the lattice spacing, a = 0.0883 ± 0.0006 fm, and the 
renormalized up/down quark mass, mP^y^' = 3.79 ±0.08 ±0.15 MeV, are consistent 
within the errors with the findings of Refs. [23j and [50j obtained at the same value 
of (3 (= 3.9). The values 2B = 4.90±0.09±0.20 GeV and F = 122.5±1.0±1.0 MeV 
correspond to a light-quark condensate equal to 

(gg)^(2 GeV) = (-264 ± 3 ± 5 MeV)=' ; (96) 

moreover, the ratio fP^y^ /F is equal to 

fphy^-/F = 1.067 ± 0.009 ± 0.009 . (97) 

The findings ( l96l) and ( 1971) are in agreement with the corresponding values obtained 
by the scaling analysis of Ref . [3l] . 

Using for the LEC's the values given in Table [9] we have calculated the values 
of the pion form factor at the physical point for the various values of considered 
in this work. Our results, including both statistical and systematic uncertainties, 
are collected in Table [10] and shown in Fig. [23], where they are also compared with 
available experimental data from Refs. [H [31 [1], [5l [6] [7] [H] • 
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(GeV") 


Fphys.fQ2\ 


0.050 


914 ± 005 ± 003 


0.148 


0.774 ±0.013 ±0.008 


0.299 


0.618 ±0.019 ±0.013 


0.503 


0.487 ±0.022 ±0.017 


0.794 


0.437 ±0.030 ±0.026 



Table 10: Values of the pion form factor F^^^^^Q"^) , extrapolated to the physical 
point using for the LEC's the results of Table\^ for various values ofQ"^. The first 
error is statistical and the second one systematic. 
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• ETMC (ChPT fit @ NNLO) 
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_L 
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0.2 
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(GeV^) 



0.6 



0.8 



Fig. 23: Pion form factor Ft, {Q'^) versus = —q^ in physical units. The full dots 
are the NNLO ChPT results of Table\T^ obtained at the physical point using for 
the LEC's the values given in Table [PI The uncertainties of the ETMC results, 
illustrated also by the dashed lines, represent the statistical and the systematic 
errors of TableHU added in quadrature. The open dots, squares and diamonds are 
experimental data from Refs. Jj^, J3, \^ and |3, 0> 13^; respectively. 



It can be seen that our values are fully consistent with the experimental data in 
the whole range of values of considered in this study. The agreement is partic- 
ularly remarkable at low values of (Q^ < 0.15 GeV^), where the experimental 
data are very precise, as well as at larger values of {Q"^ > 0.3 GeV^), where the 
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uncertainties of our results become competitive with the experimental errors. 

8 Conclusions 

We have presented a lattice calculation of the electromagnetic form factor of the 
pion obtained using the tree-level Symanzik improved gauge action with two flavors 
of dynamical twisted Wilson quarks. 

The simulated pion masses range from ~ 260 to ~ 580 MeV and the lattice 
box sizes are chosen in order to guarantee that M^L > 4. 

Accurate results for the form factor are obtained using all-to-all quark propa- 
gators evaluated with the stochastic procedure of Ref. [30] . 

The momentum dependence of the pion form factor is investigated up to values 
of the squared four-momentum transfer — 0.8 GeV^ and, thanks to the use of 
twisted boundary conditions, down to ^ 0.05 GeV^. The Q^-dependence at 
the simulated pion masses is well reproduced by a single monopole ansatz with a 
pole mass lighter by ~ 10% -h 15% than the lightest vector meson mass. 

Volume and discretization effects on the form factor have been directly eval- 
uated performing simulations at different volumes and lattice spacings, and they 
turn out to be within the statistical errors. A more complete investigation of the 
scaling properties of the pion form factor, based on the study of its mass depen- 
dence at two additional values of the lattice spacing is however desirable. The 
corresponding measurements are in progress. 

The extrapolation of our results for the pion mass, decay constant and form fac- 
tor to the physical point has been carried out using (continuum) ChPT at NNLO 
[llj. The extrapolated value of the (squared) pion charge radius is (r^)^^^* = 
0.456 ± O.OSOg^g^^ ± 0.024gyg^ in nice agreement with the experimental result 
lyiY^v- = 0.452 ± 0.011 fm^ [22j. The extrapolated values of the pion form factor 
agree very well with the experimental data up to ~ 0.8 GeV^ within uncertain- 
ties which become competitive with the experimental errors for > 0.3 GeV^. 

The relevant low-energy constants appearing in the chiral expansion of the pion 
form factor are extracted from our lattice data adding only the experimental value 
of the pion scalar radius [11] in the fitting procedure. We get: ii = — 0.4±1.3±0.6, 
I2 = 4.3 ±1.1 ±0.4, I3 = 3.2 ±0.8 ±0.2, £4 = 4.4 ±0.2 ±0.1, 4 = 14.9 ± 1.2 ± 0.7, 
where the first error is statistical and the second one systematic. Our findings are 
in nice agreement with the results of the NNLO ChPT analysis of vr — tt scattering 
data of Ref. [33]. The values found for the LEC's £3 and £4 are consistent with the 
corresponding results of the ETMC analysis of Ref. . This is quite reassuring 
because different kinds of systematic uncertainties may affect the two analyses: the 
present one being a NNLO analysis limited mainly to data from a single lattice 
spacing, and that of Ref. \3^ having two values of the lattice spacing but limited 
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mainly to a NLO analysis. 

It is the aim of our collaboration to reduce as much as possible all the un- 
certainties of the extracted low-energy constants in the next future. To this end, 
data at more values of the lattice spacing and calculations of other physical quan- 
tities, like e.g. the pion scattering lengths, will be considered. This may allow to 
avoid any input from experiments obtaining a first principle computation of the 
low-energy constants. 

In this respect a very interesting strategy is to include lattice data for the scalar 
form factor of the pion, because almost the same low-energy constants enter the 
chiral expansion of both vector and scalar form factors [TT]. In this way the use 
of the experimental value of the pion scalar radius in the fitting procedure can 
be avoided. However the lattice calculation of the scalar form factor requires the 
evaluation of both connected and disconnected diagrams. While the former have 
been already calculated on the ETMC gauge configurations, a precise evaluation 
of the latter is in progress. The results will be reported elsewhere. 
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Appendix 

In this Appendix we provide the values of the pion form factor Ft^IQ"^) obtained 
for all the simulations (see Table [T]) and for the various values of the squared 
four-momentum transfer = —q^ considered in this work. 
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